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ABSTRACT 
This  work  is  a  continuation  of  plasma  oscillation  studies 
made  in  rare  gas  plasmas  extending  over  the  last  four  years  at  the 
Naval  Postgraduate  School.   An  electronic  method  of  obtaining  con- 
tinuous graphical  plots  of  probe  current  vs.  applied  probe  potent- 
ial and  probe  current  squared  vs.  applied  probe  potential  was 
developed  and  successfully  used.   Preliminary  experiments  indicate 
this  system  may  be  extended  to  obtain  the  first  and  second  deriva- 
tives of  probe  current  with  respect  to  applied  probe  potential 
allowing  an  accurate  determination  of  space  potential  and- electron 
energy  distribution  functions.   A  gas  diode  filled  with  argon  at 
S.9  microns  of  mercury  was  used  in  this  investigation.   Probe  data 
was  taken  for  twenty  positions  along  the  axial  length  of  the  tube 
using  an  anode-cathode  spacing  of  M-2  mm.   Five  different  discharge 
currents  were  studied.   A  second  run  was  made  for  ten  positions 
and  a  discharge  current  of  50  ma  at  6.9  microns  of  mercury  to  study 
distribution  functions.   Electron  oscillations  were  detected  which 
corresponded  with  the  Langmuir  plasma  frequency  predicted  from  the 
plasma  density;  the  appearance  of  these  oscillations  corresponded 
with  a  decreasing  electron  density  gradient.   Peak  intensity  of 
oscillations  occurred  along  a  decreasing  density  gradient. 
Oscillations  were  also  measured  in  the  kilocycle  range  which  are 
attributed  to  ion  oscillations  associated  with  beam  electrons. 
Data  is  presented  in  the  form  of:   raw  data  plots;  space  potential 
as  a  function  of  probe  position;  electron  temperature  as  a  function 
of  probe  position;  electron  density  as  a  function  of  probe  position; 
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electron  density  as  a  function  of  discharge  current;  and  electron 
energy  distribution  functions.   It  is  tentatively  suggested  that 
the  unmodulated  beam  of  charged  particles,  coming  from  the  cathode, 
passing  through  and  maintaining  the  plasma,  suddenly  becomes  un- 
stable and  gives  rise  to  electron  oscillations  of  the  plasma  in 
the  presence  of  an  associated  decreasing  plasma  density  gradient. 
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1.   Introduction. 

The  recent  successes  in  harnessing  the  explosive  fusion  of 
the  nuclei  of  the  light  elements  and  the  even  more  recent  attempts 
at  space  exploration  has  led  to  a  renewed  peak  of  interest  in 
plasma  physics  in  an  effort  to  better  understand  the  basic  physical 
processes  occurring  in  a  highly-ionized,  gaseous,  conducting  medium, 

Success  in  understanding  these  phenomena  may  ultimately  lead 

to  generation  of  fusion  power  through  the  medium  of  a  plasma. 

Linhart  (1)  describes  this  possibility  in  the  following  way. 

It  was  appreciated  that  the  energy  output  from  fusion 
reactions  depends  critically  on  the  kinetic  energy  of 
the  colliding  nuclei  and  that  fusion  outputs  of 
practical  interest  depend  upon  one's  ability  to  produce 
temperatures  of  at  least  several  million  degrees  Kelvin. 
If  explosions  are  to  be  avoided,  then  the  pressure  of 
matter  at  this  temperature  must  be  balanced  by  external 
forces.   This  is  within  the  power  of  our  engineers  only 
if  the  density  of  the  nuclear  fuel  is  less  than  the 
density  of  our  atmosphere.   The  search  for  a  mechanism 
of  a  controlled  release  of  fusion  energy  became,  there- 
fore, synonymous  with  the  study  of  high  temperature, 
low  density  plasmas. 

Since  99.9%  of  the  matter  in  our  universe  is  ionized  and 

therefore  composed  of  plasma,  communication  and  electronic  systems 

used  in  space  exploration  must  be  able  to  operate  in,  and  propagate 

signals  through,  ionized  regions  of  the  outer  atmosphere  and  space . 

Thus,  an  understanding  of  plasma  processes  will  contribute  strongly 

to  the  solution  of  these  space  problems. 


Plasma  physics  is  concerned  with  the  behavior  of  systems  of  many 
free  electrons  and  ionized  atoms  where  the  mutual  Coulomb  interactions 
cannot  be  ignored.   In  a  restricted  sense,  such  systems  of  particles 
consist  of  nearly  equal  numbers  of  positive  ions  and  negative 
electrons.   Systems  of  this  type  are  examples  of  a  medium  known  as 
plasma  which  in  many  respects  behaves  differently  from  the  solid, 
liquid  and  gaseous  states  of  matter  (1)  . 

The  average  kinetic  energy  per  particle  determines  the  state  in 
which  matter  exists.   Matter  can  exist  as  a  plasma,  i.e.,  in  its 
fourth  state,  when  the  kinetic  energy  per  plasma  particle  exceeds 
the  ionizing  potential  of  atoms  --  which  is  usually  a  few  electron 
volts  (1)  . 

Plasma  spans  a  broader  energy  band  than  any  other  state  of 

matter;  it  occupies  about  20  octaves  of  the  kinetic  energy  spectrum- 

This  wide  kinetic  energy  band  of  the  plasma  state  is  the  reason  for 

so  much  common  ground  between  plasma  physics  and  many  other  fields  of 

physics,  such  as  the  dynamics  of  single  charged  particles,  or  the 

physics  of  electrical  discharges  in  gases  (1) .   Linhart  says: 

Although  probably  more  than  99.9%  of  the  matter  in  our  universe 
is  ionized  and  therefore  in  the  plasma  state,  on  our  planet 
plasma  has  to  be  generated  by  special  physical  processes  and 
special  conditions. 

On  earth  the  medium  which  most  closely  resembles  an  ideal  plasma 

is  a  partially  ionized  gas .   This  fourth  state  of  matter  defined  as  a 

plasma  is  not  very  much  at  home  on  our  dense,  cold  planet;  therefore, 

special  conditions  must  be  realized  in  order  to  generate  plasma-like 

conditions  in  the  laboratory.   These  special  conditions  are  in  the 
realm  of  the  physics  of  electrical  discharges  in  gases  (gaseous 
electronics)  and  hence,  the  close  kinship  between  this  discipline  and 
plasma  physics. 


Plasma  was  originally  defined  by  Langmuir  and  Tonks  (2) 

in  the  following  way: 

The  word  'plasma'  will  be  used  to  designate  that  portion 
of  an  arc  type  discharge  in  which  the  densities  of  ions 
and  electrons  are  high  but  substantially  equal.   It 
embraces  the  whole  space  not  occupied  by  sheath. 

The  positive  columns  of  the  glow  and  arc  discharge  are  ex- 
amples of  partially  ionized  plasmas.   Plasma  regions  are  highly 
conducting.   In  rare  gas  plasmas,  the  negative  carriers  are 
electrons.   Ionization  of  the  plasma  is  produced  essentially  by 
electron  collisions  with  some  photoionization  also  occurring. 
The  positive  ions,  the  electrons  and  the  neutral  gas  atoms  may  or 
may  not  be  in  thermal  equilibrium.   The  temperature  of  the  positive 
ions  is  usually  greater  than  the  gas  temperature,  while  the  elect- 
ron temperature  may  be  very  high.   Generally,  positive  ions  and 
electrons  may  be  assumed  to  have  a  Maxwellian  distribution; 
however,  variations  from  this  distribution  have  been  noticed 
(3,  1,  5,  6). 

Even  though  a  plasma,  with  its  ions  and  electrons  moving 
at  random  according  to  their  characteristic  temperatures,  is 
essentially  neutral  when  considering  a  large  volume,  there  must 
be  high  fields  existing  at  points  within  the  volume.   Thus,  at 
any  point  the  field  will  vary  with  time  depending  upon  the 
instantaneous  random  configurations  of  charged  particles  surround- 
ing the  point.   Fields  of  this  type  surrounding  a  point  in  a  plasma 
may  be  called  microfields  and  because  of  them,  electrons  may  gain 
considerable  amounts  of  energy  in  passing  through  a  series  of 
accelerating  microfields. 


Another  property  of  an  ionized  gas  which  owes  its  origin  to 
space  charge  phenomena  and  high  effective  conductivity  is  that  of 
electrostatic  plasma  oscillations.   These  longitudinal  oscillation-, 
of  the  plasma  are  associated  with  the  fact  that  a  small  moment, 
displacement  of  charges  of  one  sign  with  respect  to  the  other  in  the 
plasma  will  establish  an  electric  field  (microfield)  in  such  a 
direction  as  to  oppose  the  displacement.   Inertial  effects  of  the 
displaced  charges  then  provide  a  mechanism  for  periodic  oscillat i 
about  the  point  of  zero  average  charge  displacement.   The  frequency 
of  these  oscillations  is  determined  by  the  mass  of  the  displaced 
charges  and  the  'spring  constant1  of  the  electric  field  produced  by 
the  displacement   (7) . 

A  simple  derivation  of  longitudinal  plasma  oscillations  might 
be  made  as  follows : 

We  write  Newton's  second  law  of  motion, 

(1)  d/dt(m  dr/dt)    =  F  =-eE. 

Now  considering  only  wave  propagation  longitudinally, 
i.e.,  in  the  X  direction,  this  reduces  to 

(2)  mX  +  eEx  =  0 

We  also  write  Maxwell's  Curl  equation 

(3)  VxH=    1  IT  J  +  -ff   . 

a  l 

Since  no  external  magnetic  field  is  present,  H  =  0 

Now  reducing  the  equation  to  its  scalar  component  in  the 

X  direction,  and  denoting  partial  differentiation  with  respect 

to  time  by  •  ,  we  have 

O)    0  =  -  i+  IT  neX  +  Ex  . 


Then,  integrating  and  setting  the  constant  of  integration 
equal  to  zero,  since  we  have  no  external  electric  field  present, 
we  obtain 

(5)   Ex  =   M-  TT  neX 

Substituting  (5)  into  (2) ,  we  obtain 

(G)   X  +   M-  n   ne2   X  =  Q 

m 

This  equation  gives  for  the  frequency  of  oscillation 

=   8980  n^  cps. 

Where: 

r  is  the  vector  displacement  of  the  electron  from 

the  point  of  zero  charge  displacement 

e  is  the  charge  of  the  electron 

m  is  the  mass  of  the  electron 

F  is  the  restoring  force 

E  is  the  electric  field  intensity 

H  is  magnetic  field  intensity 

J  is  the  displacement  current  density 

t  is  time 

•  is  the  partial  derivative  with  respect  to  time. 

Plasma  oscillations  may  be  either  electronic  or  ionic  in 
nature,  and  the  frequencies  are  determined  by  the  respective 
electron  and  ion  masses,  hence,  ion  oscillations  will  be  much 
lower  in  frequency  than  electron  oscillations  because  of  their  much 
greater  mass   (7)  . 

Plasma  electron  oscillations  may  transmit  energy  by  a  growing 
wave  mechanism.   A  growing  wave  results  from  the  excitation  of 
plasma  electron  oscillations  and  the  interaction  of  these  oscillations 
with  the  primary  beam.   As  the  wave  moves  longitudinally  through  the 
plasma,  it  tends  to  accumulate  the  slowly  drifting  groups  of  oscill- 
ating electrons  in  advance  of  the  wave  and  disperse  the  electron 
groups  after  passage  of  the  wave.   As  this  redistribution  takes  place, 
an  exchange  of  energy  between  electrons  and  wave  in  a  particular  phase 


relationship  increases  the  energy  in  the  wave.   The  wave  amp  I  i t ude 
grows  exponentially  reaching  the  maximum  rate  of  growth  when  plasma 
frequency  approaches  the  excitation  frequency  of  the  wave  (8) . 

Energy  of  plasma  electron  oscillations  may  also  be  transmitted 
by  fast-moving  electrons  traversing  the  oscillating  region  and  having 
their  velocities  modulated  periodically.   This  type  of  energy  transfer 
obviously  depends  upon  the  relative  velocities  of  the  oscillating  and 
non-oscillating  groups  of  electrons   (3)  . 

Experiments  in  plasma  physics  concerned  with  probe  studies  and 
plasma  oscillations  are  usually  carried  out  in  gas-filled  diodes  or 
triodes  utilizing  discharges  through  rare  gases  or  mercury  vapor. 
A  diode  filled  with  Argon  was  used  in  this  study. 

Typically,  the  voltage  profile  across  such  a  tube  from  cathode 
to  anode  consists  of  a  steep  rise  in  potential  near  the  cathode, 
followed  by  a  slow  gradual  rise  across  the  positive  column  of  the 
discharge  to  the  anode  with  a  sharp,  small  rise  or  fall  in  the 
region  near  the  anode.   The  electrons  traversing  the  interelectrode 
distance  gain  their  major  potential  rise  in  region  adjacent  to  the 
cathode  called  the  cathode  TTfalln   (9).  The  cathode  is  surrounded 
by  a  sheath  of  positive  ions  because  of  the  greater  mobility  of  the 
electrons  compared  with  the  ions.   Electrons  are  accelerated  as  they 
move  away  from  the  cathode  sheath  producing  more  ionizations  and 
attaining  a  Maxwellian  distribution  of  velocities  as  they  travel 
through  the  gas  making  collisions  with  gas  atoms  and  each  other. 
At  the  negative  glow  boundary,  the  entering  electrons  are  composed 
of  two  or  more  groups  --  fast  electrons  which  have  been  produced  at 
the  cathode  and  have  experienced  few  collisions  in  the  cathode  dark 


space,  and  slow  electrons  which  have  made  many  inelastic  collisions 
in  the  dark  space.   The  slow  electrons  have  energies  below  the 
ionization  potential  of  the  gas,  but  at,  or  above,  the  excitation 
potential;  therefore,  they  make  exciting  collisions  with  the  gas 
molecules,  causing  light  emission  in  the  negative  low  region  as  the 
excited  atoms  emit  quanta  in  returning  to  their  ground  states. 
Afterward,  these  slow  electrons  have  such  low  energies  that  they 
then  recombine  with  positive  ions.   The  group  of  fast  electrons  are 
charge  carriers  which  enter  into  the  positive  column  part  of  the 
discharge.   Some  of  these  electrons  may  carry  through  to  the  anode 
depending  upon  the  mean  free  path,  but  more  likely,  these  electrons 
interact  in  the  positive  column,  give  up  some  of  their  energy  and 
become  randomized.   The  anode  will  then  receive  principally  elect- 
rons coming  from  the  positive  column. 

This  investigation  was  originally  started  at  the  United  States 
Naval  Postgraduate  School  in  1958  by  Webb  and  Garner,  continued  by 
Flynn  and  Shoults  in  1959,  and  by  Culwell  and  Johnson  in  1960.   The 
present  work  is  principally  a  continuation  of  the  work  of  Culwell 
and  Johnson  who  redesigned  and  rebuilt  the  vacuum  system  and  dis- 
charge tube,  but  were  able  to  carry  out  only  limited  rudimentary 
measurements  because  of  the  press  of  time. 

The  principal  effort  of  this  work  was  in  (a)  the  design  and 
development  of  a  suitable  electronic  system  with  which  to  make 
intensive  probe  studies  in  the  discharge  tube  and  (b)  then  the 
making  of  such  measurements.   A  unique  system  was  designed  yielding 
continuous  graphical  plots  of  conventional  Langmuir  data  such 


probe  current  vs.  applied  probe  potential  and  probe  current 
squared  vs.  applied  probe  potential.   This  system  shows  promise 
of  extension  to  yield  accurate  space  potential  determination  and 
electron  energy  distribution  functions. 

This  report  covers  plasma  oscillation  and  probe  studies  made 
in  low  pressure  argon  discharges. 


2 .   Previous  Work 

Investigation  of  gaseous  discharge  phenomena  at  low  pressures 
was  given  its  first  real  impetus  in  192M-  when  Langmuir  and  Mott- 
Smith  (10)  reported  a  new  method  of  making  extensive  studies  of 
electrical  discharges  in  gases  at  low  pressures.   Their  technique 
consisted  of  placing  a  small  auxiliary  electrode  (commonly  referred 
to  as  a  probe)  in  the  discharge   and  studying  the  complete  probe 
volt-ampere  characteristic  using  the  collector  theory  proposed 
by  Langmuir  (11)  . 

Further  studies  by  Langmuir  (12)  in  1925  disclosed  that  appar- 
ently there  were  three  energy  groups  of  electrons  present  in  the 
discharge.   He  classified  these  groups  as  "primary  electrons" 
(those  electrons  which  keep  nearly  all  the  momentum  gained  in  the 
cathode  ion  sheath) ,  "secondary  electrons"  (those  electrons  having 
a  Maxwellian  distribution  of  energies  proportional  to  the  energy  of 
the  primaries) ,  and  a  third  group  of  "ultimate  electrons"  (those 
electrons  having  a  Maxwellian  distribution  of  energies  independent 
of  the  primaries  and  much  lower  than  the  secondaries) . 

These  observations  indicated  that  there  must  be  some  type  of 
energy  transfer  mechanism  involved  to  permit  certain  electrons  to 
gain  energies  in  excess  of  the  tube  drop  while  others  became  energy 
deficient.   Langmuir  referred  to  this  phenomenon  as  scattering. 


Floating  probes  had  been  used  previously  in  analyzing  discharges, 
however,  they  were  of  limited  value  because  of  the  failure  to  appre- 
ciate sheath  formation  surrounding  the  probe  and  the  consequent 
failure  to  measure  space  potential. 


Dittmer  (13)  experimentally  disproved  the  hypothesis  that  a 
new  mechanism  (such  as  the  possible  large  increase  in  the  target 
area  of  excited  molecules)  was  the  cause  of  scattering  by  showing 
that  the  scattering  reached  a  constant  value  at  some  distance  from 
the  cathode.   Dittmer  attempted  to  find  oscillations  in  an  effort 
to  validate  the  hypothesis  that  electrical  oscillations  in  the  dis- 
charge caused  scattering.   He  was  unsuccessful  in  all  attempts  at 
finding  such  oscillations. 

Tonks  and  Langmuir  (2)  developed  a  simple  theory  for  electron 
and  ion  oscillations  which  was  analogous  to  the  electrons  being  im- 
bedded in  a  jelly-like  substance  composed  of  ions.   The  electrons 
have  a  natural  frequency  given  by 


fP   =  /ne2  ]      cps 
[  tf  me/ 

where  n  is  the  electron  density  per  cubic  centimeter,  e  the  electronic 
charge,  me  the  mass  of  the  electron.  The  ions  obey  a  similar  law  with 
a  natural  frequency  approaching 

f±      =   fe   (me/nO 
where  m-  is  the  mass  of  the  ion. 

Tentative  verification  of  this  theory  was  obtained  by  detecting 
oscillations  in  the  300  to  1000  megacycle  range  which  corresponded 
reasonably  well  with  the  predicted  frequencies.   It  was  conjectured 
that  additional  oscillations  in  the  50  to  60  megacycle,  1.5  megacycle 
and  lower  ranges  might  correspond  to  oscillations  of  the  primary 
beam  electrons  and  ions,  respectively. 
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Merrill  and  Webb  (1M-)  carried  on  extensive  investigat  Lons  in 
a  mercury  arc  discharge  at  low  pressures  and  discharge  currents 
with  a  movable  probe  to  determine  a  possible  relation  between  the 
plasma  oscillations  and  scattering  as  suggested  by  Tonks  and 
Langmuir.   They  found  that  scattering  occurred  in  several  narrow 
regions  and  was  associated  with  periodic  oscillations  of  consid- 
erable magnitude  which  agreed  with  the  Langmuir  frequency  within 
ten  per  cent.   The  fact  that  several  regions  were  observed  has 
been  explained  by  Bohm  and  Gross  to  be  a  bunching  effect  of  the 
electron  beam.   Merrill  and  Webb  further  noted  that  for  each 
scattering  region  there  was  only  one  observable  frequency  and 
that  the  maximum  intensity  of  the  oscillations  occurred  at 
from  one  half  to  one  millimeter  from  the  scattering  region  on 
the  anode  side. 

Merrill  and  Webb  concluded  that  the  scattering  was  due  to 
the  oscillations  receiving  energy  from  the  fast  electrons,  which 
then  give  some  of  the  electrons  excess  energies,  and  that  the 

observed  difference  in  the  positions  of  the  oscillations  and  the 

* 
scattering  was  a  consequence  of  the  experimental  technique. 


•k 

With  the  recent  theory  proposed  by  Allis  and  partially 
confirmed  by  Emeleus,  one  might  find  that  an  analysis  of  Merrill 
and  Webb's  600  complete  volt-ampere  characteristic  curves  would 
show  that  their  observation  between  scattering  and  oscillation 
locations  was  of  considerable  physical  significance. 
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Although  no  major  breakthrough  in  understanding  oscillations 
occurred  in  the  next  ten  years,  many  important  general  observation^ 
were  made.   A  few  of  these  are  stated  below  (15) . 

(a)  The  operating  position  on  the  current  -  voltage  char- 
acteristic of  the  tube  under  investigation  may  affect 
the  amplitude  of  oscillations  observable  in  external 
circuits.   It  has  been  found  that  oscillations  can 
usually  be  detected  over  wide  sections  of  the 
characteristic  for  oxide-coated  cathodes. 

(b)  The  position  of  the  probe  in  the  discharge  markedly 
affects  the  amplitude  of  the  observed  oscillations; 
however,  the  observed  frequency  usually  remains 
about  the  same. 

(c)  The  signal  amplitude  of  the  oscillations  depends 
considerably  on  the  electrode  geometry  and  upon 
which  electrodes  are  used  for  detection  of 
oscillations . 

Emeleus  (16)  has  pointed  out  the  value  of  visual  observations 
in  electrical  discharge  studies  using  oxide-coated  cathodes  in  ar- 
gon at  a  few  microns  of  pressure  with  a  low  voltage  discharge.   A 
representative  sketch  of  a  gaseous  discharge  in  which  high  frequen- 
cy oscillations  are  present  may  be  shown  as  follows: 
12       3         4   _     A 
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Fig.  A.   C,  cathode;  1,  cathode  sheath;  2,  unscattered  beam 
section  of  discharge;  3,  meniscus;  M-,5,  scattered  beams;  A,  anode 
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The  principal  static  fall  in  potential  occurs  in  region  1; 
the  rest  of  the  discharge  is  essentially  at  anode  potential. 
Between  1  and  2  represents  the  familiar  Dittmer  scattering  region 
with  dispersion  occurring  in  or  just  after  3.   It  has  further  been 
noted  that  strong  oscillations  have  only  been  detected  to  the 
right  of  3 . 

Prior  to  1950  little  attention  had  been  given  to  the  effect 
of  probe  contamination  on  the  validity  of  the  volt-ampere  probe 
characteristic.   Easley  (17)  in  probe  studies  of  a  mercury  dis- 
charge buffered  with  argon  showed,  by  drawing  sufficient  electron 
current  to  make  the  probe  cherry- red  (and  thus  outgas  it  by  ohmic 
heating)  before  each  reading  that  many  of  the  non-linear  logarith- 
mic plots  previously  obtained  were  due  to  a  change  in  contact 
potential  and  thus  did  not  indicate  the  presence  of  a  non-Maxwell- 
ian  distribution  of  electron  velocities,  but  rather  probe  contam- 
ination. 

The  results  of  Oleson  and  Found  (18)  have  shown  that  non- 
linear logarithmic  plots  were  obtained  in  discharges  in  which 
moving  striations  were  present.   The  logarithmic  plots  obtained 
looked  as  though  two  Maxwellian  distributions  were  present  in 
the  discharge.   However,  it  had  become  fairly  evident  that  in 
many  cases  (such  as  in  discharges  at  pressures  of  the  order  of 
several  microns  of  mercury)  that  a  Maxwellian  distribution  of 
energies  did  not  exist.   Consequently,  increased  effort  was 
directed  toward  new  techniques  of  determining  the  actual  energy 
distribution  function  rather  than  assuming  it  to  be  Maxwellian. 
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It  is  worthwhile  to  note  that  Langmuir  considered  the  energy 
distribution  function  problem  in  his  original  theory  on  collectors, 
and  though  he  did  not  employ  the  technique,  he  did  show  that  ob- 
taining the  second  derivative  of  probe  current  with  respect  to 

* 

probe  voltage  would  give  the  true  energy  distribution. 

Haigh  and  Smith  (19)  undertook  a  study  of  the  energy  distri- 
butions in  the  positive  columns  of  low  pressure  (.005  mm  Hg  + 
.085  mm  Hg) ,  hot  cathode  discharges  in  argon  and  air  by  a  magnetic 
spectrum  method.   A  continuation  of  this  study  by  Haigh  and 
McCulloch  (20)  obtained  data  on  air,  argon  and  nitrogen.   Their 
results,  using  argon,  are  shown  in  Figs.  B  through  E.   These 
plots  illustrate  the  types  of  distributions  obtained  in  argon. 


Pressure  0.0825  mm  Hg 
Tube  Current  -  15  ma 


Pressure  0.026  mm  Hg 
Tube  Current  -  15  ma 
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Fig.    B.  Fig.    C. 

Strictly  speaking,  this  is  true  for  a  cylindrical  collector  if 
the  velocity  distribution  is  isotropic.   In  cases  in  which  one 
doubts  the  presence  of  an  isotropic  distribution,  a  spherical 
collector  is  required.   This  is  shown  in  Langmuir' s  original 
development. 
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Pressure  0.009  mm  Hg 
Tube  Current  -  15  ma 


Pressure  0.003  mm  Hg 
Tube  Current  -  IS  ma 
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Fig.    D. 
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Fig.    E. 

It  is  interesting  to  note  that  for  a  given  pressure,  a  decrease 


in  current  has  the  same  effect  as  a  decrease  in  pressure  at  constant 
current.   Haigh  and  McCulloch  hypothesized  that  the  second  peak  in 
their  distributions  may  be  the  result  of  the  "combined  effects  of  the 
Ramsauer  free  paths  and  the  tendency  for  a  larger  proportion  of  high 
velocity  electrons  as  the  pressure  decreases." 

As  far  back  as  1939  Druyvesteyn  (21)  had  shown  that  the  electron 
energy  distribution  function  could  be  found  from  the  expression 


'2m   V 
N  f  (energy)  =N0(eV)  =  _2_        P 

Ae 


T 


dVr 


where  N0(eV)  is  the  number  of  particles  per  unit  of  energy  per  unit 
volume,  and  in  this  work,  customarily  put  in  units  of  reciprocal 
electron  volts  and  reciprocal  cubic  centimeters,  i.e.,  eV   cm   ; 
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A  is  the  area  of  the  probe,  m  is  the  mass  of  the  electron, 
e  is  the  electronic  charge,  i   is  the  probe  current  and  V  is 
the  probe  voltage.   All  are  expressed  in  consistent  units. 

Despite  this,  no  significant  effort  was  expended  to  obtain 
energy  distributions  because  of  the  inherent  inaccuracies  in- 
volved in  the  tedious  graphical  method  of  obtaining  the  second 
derivatives  of  probe  current  with  respect  to  probe  voltage. 

Sloane  and  MacGregor  (22)  developed  an  idea  originally  pro- 
posed by  Emeleus  for  electrically  obtaining  the  energy  distribution 
By  putting  a  small  alternating  signal  on  the  probe  it  was  shown 
that  the  dc  increase  caused  by  the  ac  signal  was  directly  propor- 
tional to  the  second  derivative  of  probe  current  with  respect  to 
probe  voltage.   Though  it  appeared  that  this  technique  might  yield 
the  first  real  gain  in  obtaining  distribution  functions,  Emeleus 
mentioned  in  Loeb  (3)  that  "not  too  much  confidence  should  be 
placed  in  the  distribution  derived  by  this  technique 
principally  because  of  the  many  ways  the  plasma  might  react  to  the 
alternating  potential. 

Recently,  Boyd  and  Twiddy  (6)  used  an  improved  version  of  the 
Sloane  and  Mac  Gregor  technique  by  amplitude  modulating  that  ac 
signal  thus  permitting  work  in  plasma  having  higher  noise  levels. 
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A  typical  energy  distribution  measurement  by  Boyd  and  Tw.iddy 
for  a  striated  head  in  hydrogen  shows  at  least  two  electron  energy 
groups,  as  indicated  in  Fig.  F. 
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Fig.  F. 
Boyd  and  Twiddy  have  also  shown  that  for  a  non-isotropic, 
strongly  directed  discharge  a  cylindrical  probe  will  give  an  energy 
distribution  which  will  have  negative  values  of  energy  between  two 
major  electron  groups.   In  every  case  studied  by  Boyd  and  Twiddy, 
two  well-separated  groups  of  electrons  (and  sometimes  a  third  group) 
were  detected.   Further,  it  has  been  suggested  that  the  energy  ex- 
change between  the  electron  groups  might  be  a  much  more  rare  phenom- 
enon than  had  been  expected.      -,y 


In  1956,  Medicus  (23)  developed  a  simple  technique  for  deter- 
mining the  energy  distribution  function  from  the  probe  character- 
istic curves  plotted  automatically  with  an  X  -  Y  recorder. 
Medicus'  method  is  a  graphical  one  in  which  tangents  are  drawn  to 
the  probe  curves  and  the  differences  between  the  tangents  repre- 
sent the  second  derivatives.   Fig.  G  below  outlines  the  method 
pictorially , 
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Fig.  G 


To  obtain  the  best  results  the  probe  curve  should  be  ampli- 
fied so  that  all  regions  of  interest  lie  at  about  a  45°  slope. 
The  actual  probe  curves  may  be  used  directly  only  if  the  ion  current 
is  small  compared  to  the  electron  current,  or  if  the  ion  current 
varies  linearly  with  voltage;  otherwise,  one  must  use  the  net 
electron  current. 
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Medicus  (24)  obtained  results  which  seem  to  further  substanti- 
ate the  suggestion  of  Boyd  and  Twiddy  concerning  the  possible  low 
energy  exchange  between  slow  and  fast  energy  groups.   Medicus  showed 
that  when  a  discharge  consists  of  both  a  high  energy,  non-Maxwellian 
group  of  electrons  and  a  low  energy  Maxwellian  group  that  "the 
characteristic  linear  part  of  the  corresponding  probe  curves,  extra- 
polated to  plasma  potential,  directly  determines  the  random  current 
density  of  the  fast  group."   The  expression  for  this  relationship 


is  given  by   4-  iT 
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Fig.  H  shows  how  this  is  obtained  graphically  and  how  it  is  re- 
lated to  the  high  energy  electron  distribution. 
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Applying  the  diffusion  law  to  his  result,  Medicus  found  th.it 
the  interaction  between  the  slow  and  fast  groups  must  be  small  and 
that  cases  may  exist  "at  low  total  plasma  densities  where  the 

particle  density  of  the  slow  group  is  always  comparable  to,  or 

* 
even  much  higher  than,  that  of  the  fast  group." 

Continuing  the  search  to  improve  probe  techniques  considerable 
effort  has  been  expended  to  develop  efficient  methods  of  eliminating 
probe  contamination  while  obtaining  the  probe  characteristic  curves. 

Wehner  and  Medicus  (25) ,  using  an  automatic  X  -  Y  plotter, 
were  able  to  obtain  a  complete  probe  characteristic  curve  in  sever- 
al seconds  coming  from  a  cherry-red  probe  condition.   They  claimed 
no  deviation  from  a  Maxwellian  distribution  over  nearly  M-  powers 
of  10  in  probe  current  even  in  tubes  with  oxide-coated  cathodes, 
providing  the  characteristic  was  recorded  in  several  seconds  after 
the  probe  had  been  thoroughly  cleaned  as  described  above. 

White  (26) ,  while  making  measurements  in  a  hollow-cathode 
glow  discharge,  found  that  excessive  material  was  sputtering  from 
the  cathode  onto  the  probe.   It  was  found,  that  by  keeping  the 
probe  at  30  volts  below  floating  potential  and  by  pulsing  it 
every  0.1  seconds  for  0.004  seconds,  it  was  possible  to  maintain 
a  clean  probe  by  ion  bombardment  and  hence,  obtain  a  reliable 
probe  characteristic. 


This  would  imply  that  the  particle  densities  of  the  different 
energy  groups  would  have  discrete  values  at  any  point  in  the 
discharge  tube  for  a  given  set  of  parameters. 
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Several  months  later  Waymouth  (27)  reported  d  similar 
pulse  technique  for  obtaining  "uncontaminated"  probi 
characteristics . 

Some  results  of  electron  current  vs.  probe  voltage  for 
discharges  at  pressures  of  six  microns  of  mercury  plus  one 
mm  of  neon  are  shown  in  Fig.  I,  below,  taken  from  Waymouth ?s 
article  (27) . 
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An  attempt  to  obtain  energy  distribution  functions  by  the 

method  of  Medicus  gave  poor  results  because  of  the  small  size  of 

the  oscilloscope  photographs.   Standard  scatter  of  data  is  shown 

for  a  probe  0.25  cm  from  the  cathode. 
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It  was  found  that  this  pulse  technique   represents  at 
least  a  600-fold  reduction  in  the  net  rate  of  probe  cont.im 
ination  in  gas  discharges.   It  was  also  shown  that  the  energy 
distribution  of  electrons  in  the  negative  glow  consists  of 
two  groups,  a  Maxwellian  group  and  a  high  energy  group  with  a 
potential  approximately  equal  to  the  cathode  fall. 


It  is  important  that  in  the  zeal  to  overcome  one  problem 
it  is  not  replaced  with  a  more  complex  one.   For  these  measure- 
ments to  have  meaning  the  ion  sheath  must  be  in  equilibrium; 
therefore,  it  is  clear  that  there  is  some  upper  limit  on  the 
rise  time  of  the  pulse  so  that  the  sheath  thickness  is  essenti- 
ally in  constant  equilibrium,  i.e.,  where  the  distance  the  ions 
can  travel  during  a  pulse  equals  the  product  of  the  pulse  dura- 
tion, the  ion  mobility  and  the  sheath  field. 
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Waymouth  (28)  has  pointed  out  an  important  rn 
diagnostic  tool  in  making  a  rapid  evi 
principal  electron  energy  groups  directly  from  the  linear 

i  vs  W,  probe  characteristic  curves.   His  technique  is 

p     P  ^ 

essentially  this . 


h  ,2 


N0(eV)  ~  V^  dip/dVp 


therefore,  whenever  the  linear  ip  vs  V  curve  has  a  constant 
slope 

N0(eV)  =  0, 
i.e.,  there  are  no  electrons  of  corresponding  energy  in 
the  energy  distribution  spectrum. 

An  example  is  shown  in  the  following  plots: 


E(ev) 


Space 

Potential 


Fig.  K 
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A  modern  radio  frequency  technique  of  determining  the 
plasma  resonant  frequency  in  a  gaseous  discharge  has  be< 
developed  by  Takayama,  Ihegumi ,  and  Miyozaki  (29).   They 
discovered  that  if  an  rf  signal  is  superposed  on  a  probe, 
which  is  maintained  negative  at  all  times  relative  to  space 
potential,  then  a  resonant  current  increase  appears  in  the 


circuit  at  the  Langmuir  frequency,  f  =  (   ne 


2  \1> 


cps 


The  relation  between  the  resonance  and  the  frequency  of 
imposed  rf  voltage  is  as  shown  in  Fig.  L,  below,  for  a 
homogeneous  mercury  plasma  at  0.1  mm  Hg  pressure. 
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Takayama,  et  al.,  suggest   that  this  is  a  beam-plasma 
type  of  interaction  and  that  the  resonance  effect  will  not 
be  discernible  when  the  energy  supplied  by  the  rf  field  is 
comparable  to  the  thermal  energy  of  the  electrons.   They 
point  out  that  an  increase  in  pressure  will  likewise  cause  a 
decrease  in  resonance  since  the  height  of  the  resonance  peak 

also  depends  upon  the  collision  frequency  of  the  electrons  with 

atoms . 
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In  19S7,  Allis  (30)  suggested  that  as  an  electron  beam 
travels  through  a  region  of  increasing  plasma  density,  electron 
oscillations  are  being  generated  in  a  region  in  which  they  can- 
not propagate  in  the  absence  of  the  beam.   On  the  other  hand, 
an  electron  beam  travelling  through  a  plasma  of  decreasing 
density  might  generate  oscillations  which  could  propagate  and 
outrun  the  beam. 

Allis  further  pointed  out  that  if  a  monoenergetic  beam  of 
electrons  is  going  through  a  plasma  then  the  average  beam  velocity 
is  reduced  by  the  growing  amplitude  of  the  oscillations,  the 
energy  being  transferred  into  field  energy  in  the  wave  and  in- 
creased kinetic  energy  of  the  random  plasma  electron  group. 
Allis  states,  "growth  occurs  mainly  below  the  plasma  frequency" 
so  that  the  free  plasma  cannot  transmit  these  waves,  (see  Fig.  M) . 
He  further  mentions  that  if  the  beam  should  peter  out  through 
scattering  or  slowing  down,  the  plasma  waves  must  decay  expon- 
entially along  with  the  beam. 


Amplification 
Factor 


Fig.  M.  Amplification  factor  for  an  electron 
beam  in  a  plasma  (30) . 
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A  recent  investigation  of  electron  beam  interaction  with  a 
plasma  was  conducted  by  a  group  of  Russian  workers  (31) .   In  their 
paper  they  point  out  that  theoretical  investigations  have  shown 
that  the  state  of  an  unmodulated  beam  of  electrons  moving  through 
a  homogeneous  plasma  is  unstable.   Density  and  velocity  fluctuations 
arising  in  the  beam  or  the  plasma,  may  excite  longtudinal  electro- 
magnetic waves  (plasma  oscillations)  of  growing  amplitude  in  the 
plasma.   They  note  that  the  initial  fluctuations  may  be  caused  by 
external  forces. 

These  workers  state  that  the  observed  maximum  excited  or 
amplified  frequency  is  determined  by  the  following  relation: 


Cjl)t 


^o  w      ~P ,  for  vq  ^>> 

'l  -      S' 


2  \h 


2 
vo 

Where, 

wo  is  the  observed  frequency  in  radians 
COp  is  the  Lnagmuir  plasma  frequency 
"Q  is  the  beam  velocity 
S  is  the  average  thermal  plasma  velocity 
They  also  point  out  that  a  large  velocity  spread  in  the  beam 
may  remove  the  instability.   If  the  regularly  moving  beam  consists 
of  ions,  the  excitation  is  less,  not  by  a  factor  me/m^,  but  only 
(me/m^)    .   Therefore,  the  interaction  mechanism  considered  is 
also  important  if  an  ion  beam  moves  through  a  plasma. 
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Plasma  instabilities  may  be  either  convective  or  absoluti 
Convective  instabilities  tend  to  move  with  the  flow  and  the 
oscillations  they  produce  amplify  or  diminish  with  time.   Ab- 
solute instabilities  produce  oscillations  which  grow  contin- 
uously with  time. 

Emeleus  and  Garscadden  (32)  have  very  recently  done  some 
work  in  exploring  plasma  oscillations  with  a  transverse  electron 
beam  observing  the  resultant  energy  spread  of  the  transverse  beam 
when  it  is  in,  or  has  traversed,  a  section  of  the  main  discharge 
where  scattering  or  oscillations  can  be  detected. 

A  typical  result,  for  a  discharge  of  1.5  microns  pressure 
in  mercury  vapor  with  a  discharge  current  of  27  ma  (8  mm  diameter 
cathode)  with  an  electrode  separation  of  M-2  mm  and  a  meniscus 
forming  10  mm  from  the  cathode,  is  shown  below. 
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Fig.  N.   Energy  spread  for  a  transverse  electron 
beam  in  a  meniscus-type  discharge. 
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The  conclusion  is  made  that  "with  the  discharge  described 
*  *  *  The  results  appear  to  be  compatible  with  considerab 1 1 
randomization  and  possible  break-up  of  the  main  beam  by  the 

oscillation"  as  the  beam  interacts  with  the  plasma  in  the 

* 

meniscus  region. 

Very  recent  investigations  by  Mahaffey,  McCullagh,  Garscadden 
and  Emeleus  (33)  are  pertinent  in  correlating  growth  or  decay  of 
oscillations  with  plasma  density  gradients. 

Their  work  was  also  done  in  a  mercury  discharge  at  approxi- 
mately one  micron  of  pressure  using  an  eight  mm  oxide-coated 
cathode  with  a  tube  current  of  27  ma  and  a  tube  voltage  of  about 
20  volts.   Two  main  sets  of  measurements  were  made  with  emphasis 
on  (a)  frequencies,  intensities,  and  energy  spread  of  the  electron 
beam  (found  from  the  break  in  the  i  vs.  Vp  curve  at  which  fa 
electrons  began  to  be  received) .   (b)  frequency  measurements, 
electron  temperatures,  and  electron  concentrations  (found  from 
in  vs  V  curves  near  space  potential) . 

They  conclude  that: 

Some  correlation  has  been  found  between  growth  or 
decay  of  oscillations  and  motion  of  the  electron 
beam  down  or  up  a  plasma  concentration  gradient, 
which  appears  to  be  in  general  accord  with  the 
plasma  dispersion  theory. 


*cf.  our  results,  Section  5  of  this  report,  for  a  beam  generated 
in  the  discharge  tube  interacting  with  the  plasma. 
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3  .   Equipment . 

The  equipment  used  can  be  divided  into  four  main  componcni 
groupings  as  follows: 

(1)  Gaseous  discharge  tube  in  which  a  glow  discharge  could 
be  continuously  maintained. 

(2)  Vacuum  system  for  evacuation  and  filling  of  the  tube . 

(3)  Frequency  detection  and  measuring  equipment. 

(•4)   dc    electrical  circuit  for  supplying  power  and  for 
controlling  and  measuring  discharge  characteristics. 

Figs.  1  through  3  are  photographs  of  the  general  experi- 
mental setup. 

a.   Discharge  Tube .   Figs.  4a  and  4b  are  photographs  and  Fig. 
5  is  a  detailed  drawing  of  the  tube  (diode) .   This  is  the  same 
tube  used  by  Culwell  and  Johnson  in  their  work  during  1960.   No 
changes  were  made  to  the  tube,   The  electrodes  and  probes  were 
attached  to  soft  iron  cylindrical  cores  sealed  in  glass  which 
made  possible  the  variation  of  interelectrode  distances  and 
probe  positions  by  moving  the  iron  cores  with  externally 
applied  magnetic  fields. 

The  probe  used  in  this  work  was  a  cylindrical  probe  made 
of  .008  inch  diameter  tungsten  wire  enclosed  in  a  glass  sheath 
up  to  3.8  millimeters  from  the  tip. 

The  ring  probe  shown  in  Fig.  5  was  not  used  during  the 
experiment.   It  was  damaged  during  electrode  movement,  breaking 
loose  from  its  mounting  and  falling  to  the  bottom  of  the  discharge 
tube  . 
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The  cathode  was  a  Philips  impregnated  "L"  cathode  with  a 
two  centimeter  diameter  face  indirectly  heated  by  a  tungsten 
element  drawing  seven  amperes  dc . 

The  anode  was  a  three  centimeter  diameter  nickel  disc. 

This  tube  was  originally  designed  with  high  purity  of  the 
gas  being  a  primary  design  criterion.   As  such,  the  coiled  wire 
leads  from  the  electrodes  and  the  probe  to  the  external  connect- 
ions were  unavoidable.   The  inductive  effects  of  these  leads  may 
be  considerable  at  electron  oscillation  frequencies  and  tend  to 
mitigate  against  making  large  amplitude  high  frequency  measure- 
ments through  them. 

b.   Vacuum  System.   The  vacuum  system  built  by  Culwell  and 
Johnson  (8)  was  not  changed.   It  remained  operative  during  the 

entire  course  of  this  work  without  leaking.   Pressures  of  the 

-8 
order  of  10"   mm  of  mercury  were  readily  attained.   All  valves 

in  the  system  were  glass  stopcocks  lubricated  with  Type  N  Apiezon 
grease.   The  diffusion  pump  was  a  three-stage  air-cooled  pump 
using  Octoil-S.   The  major  part  of  the  system  was  baked  out  with 
a  portable  oven  (see  Fig.  1  ) .   Heating  tapes  were  used  to  bake 
out  external  glass  tubing.   An  induction  heater  was  used  to  out- 
gas  the  anode.   Outgassing  of  the  Langmuir  probe  was  accomplished 
by  ohmic  heating  while  the  tube  was  operating. 


Manufactured  by  Semicon  of  California,  Inc 
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Low  pressures  were  measured  with  a  Veeco  Type  RG  31  A 
ionization  gage .   Pressures  in  the  operating  range  of  one  to 
ten  microns  were  measured  with  a  Langmuir  viscosity  gage  (35) . 
Tube  voltage  readings  were  also  used  as  an  indication  of  tube 
pressures  after  the  voltage-pressure  relation  had  been  determined, 
Tube  voltage  is  a  function  only  of  tube  pressure  if  other  para- 
meters such  as  anode-cathode  spacing  and  discharge  current  remain 
constant;  therefore,  a  given  pressure  could  be  readily  duplicated 
by  evacuating  or  filling  the  tube  until  the  correct  voltage  was 
obtained.   Thus,  when  fresh  argon  was  desired  in  the  tube,  flush- 
ing and  refilling  was  quickly  and  accurately  accomplished  without 
relying  on  the  slower  pressure  measurement  process  using  the 
viscosity  gage . 

c.      Frequency  Detection  Equipment.      A  cavity  wavemeter,    Test 
Set,    Model   AN/UPM   -    2,    with   a   frequency   range   of  80    to   1220   Mcs 
was    used   to   detect   electron  oscillation   frequencies  by  coupling 
with   the  wire  probe    at   anode  potential.      The   dc    signal   from  the 
cavity  wavemeter  detector  crystal  was   fed   into   a  Hewlett-Packard 

Model   M-25  A  dc  micro  volt-ammeter  to  obtain  detectable   signals . 

-12 
Currents  as  small  as  10    amperes  could  be  detected  with  this 

instrument . 

A  Hewlett-Packard  UHF  signal  generator,  Model  612  A,  with 

a  frequency  range  of  M-50  to  1230  Mcs  was  also  tried  utilizing  the 

technique  of  parametric  amplification  suggested  by  K.  Takayama, 

et.  al .  (29)  which  essentially  consisted  of  putting  the  rf 
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potential  on  the  probe  held  at  space  potential  while  the  freqm  ncy 
of  the  generator  was  varied  in  an  attempt  to  cause  resonance  at 
the  electron  oscillation  frequency. 

Ion  oscillations  were  detected  with  a  Tektronix  Type  581 
oscilloscope  and  also  simultaneously  counted  with  a  Hewlett-Packard 
Electronic  Counter,  Model  521  A.   The  ion  signals  were  amplified 
by  a  Hewlett-Packard  "40  db  amplifier,  Model  1+50  A,  to  feed  count- 
able signals  to  the  counter. 

d.   dc  Circuitry  Components.   Fig.  6   is  a  schematic  block  dia- 
gram showing  the  general  dc   circuitry  employed. 

Filament  voltage  was  supplied  by  storage  batteries  regulated 
by  a  charging  circuit  during  operation.   Seven  amperes  dc  was 
sufficient  to  heat  the  Philips  L  type  cathode  to  1050°  C.   Cathode 
temperature  was  measured  on  the  face  with  a  Leeds  and  Northrup 
optical  pyrometer.   Cathode  temperature  was  maintained  within  + 
5°  C  during  operation. 

Electrode  spacing  was  measured  with  a  locally  manufactured 
cathetometer  which  allowed  electrode  positioning  within  +  0.2  5  mm. 

The  discharge  tube  was  powered  by  a  Kepco  Labs  Regulated 
High  Voltage  Power  Supply,  Model  1250  B,  having  a  voltage  range 
from  0  to  1000  volts  and  a  current  output  from  0  to  500  milliam- 
peres . 

An  F.  W.  Moseley  Autograf  X  -  Y  Recorder,  Model  2,  was  used 
in  the  plate  circuit  to  obtain  continuous  plots  of  the  volt  - 
ampere  tube  characteristic  simultaneously  with  the  taking  of  probe 
data,  thus  the  effect  on  discharge  current  of  drawing  probe 
current  could  be  monitored  and  recorded. 
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In  i  LaJ      pts    taking  probe  data  were  designed  around 
oscillographic  recording  of  the  probe  characteristic  curves. 
Spurious  high  frequency  signals  in  the  range  0,5  to  5  Mcs  of  un- 
known origin  interfered  with  these  attempts  despite  coaxial 
shielding,  etc.   The  oscillographic  method  was  therefore  abandoned 
as  a  primary  means  of  measurement  and  a  system  designed  to  mini- 
mize ac  interference  was  constructed.   A  description  of  this  system 
follows . 

The  probe  voltage  was  supplied  by  a  Hewlett  -  Packard  Low 
Frequency  Function  Generator,  Model  202  A,  set  for  triangle  wave 
generation  to  obtain  a  voltage  source  varying  linearly  with  time. 
This  signal  was  then  amplified  and  biased  to  obtain  the  desired 
positive  and  negative  voltage  excursions  relative  to  the  anode . 

A  Donner  Analog  Computer,  Model  30,  with  Problem  Board,  Model 
30-6,  was  used  to  amplify  the  voltage  drop  across  the  variable 
resistor  in  the  probe  circuit,  thus,  producing  a  voltage  proportion- 
al to  probe  current  in  order  to  ensure  a  voltage  range  suitable  for 
best  operation  of  the  function  multiplier  and  X  -  Y  recorder.   Ten 
operational  amplifiers  were  available  with  this  computer.   The 
analog  computer  was  also  used  as  a  source  of  five  separate  regulat- 
ed dc  power  supplies  from  its  initial  conditions  panel.   The  analog 
computer  amplifiers  operate  in  the  range  -100  to  +100  volts  with 
long  time  stability  better  than  0.5%  full  scale  value.   Initial 
condition  voltage  sources  operate  from  0  to  100  volts.   Computer 
amplifiers  have  a  maximum  drift  characteristic  of  M-  millivolts 
per  hour,  with  a  short  term  random  drift  less  than  +  2  mv. 
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Balancing  of  amplifiers  is  easily  accomplished  and  holds  dri I  t 
error  to  negligible  levels  over  short  time  intervals  for  high 
signal  to  drift  ratios.   Potentiometers  on  the  problem  board 
allowed  precise  gain  control  on  amplified  signals. 

A  Donner  Electronic  Function  Multiplier,  Model  33,  was  used 
to  square  the  values  of  voltage  proportional  to  probe  current . 
This  equipment  would  accept  inputs  of  +  100  volts  in  either  or  both 
input  channels  (X,  Y)  and  produced  an  output  voltage  of  -0.01  XY. 

An  F.  W.  Moseley  Autograf  X  -  Y  Recorder,  Model  2S,  was  used 
to  make  continuous  plots  of  probe  data.   Probe  voltage  with  re- 
spect to  the  anode  was  applied  to  the  horizontal  axis  which  has 
a  range  from  0-7.5  millivolts  to  0-150  volts  full  scale.   The 
voltages  proportional  to  probe  current,  or  probe  current  squared, 
were  applied  to  the  vertical  axis  which  has  a  range  0-5  millivolts 
to  0  -  100  volts  full  scale.   Input  impedances  of  both  axes  are  of 
the  order  of  2  megohms.   Sweep  rate  was  controlled  by  the  low 
frequency  function  generator.   It  was  generally  kept  from  0.05 
cps  to  0.2  cps  to  avoid  inertial  error  in  the  X  -  Y  recorder. 
This  error  is  negligible  at  these  sweep  rates.   The  Autograf  is 
accurate  to  0.25%  of  full  scale  value  on  all  ranges  when  calibrated. 
Calibration  drift  is  negligible . 

Both  the  output  of  the  function  multiplier  and  input  of  the 
X  -  Y  recorder  have  low  pass  filters  keeping  ac  Interference  to 
a  minimum. 


31+ 


A  Lambda  Electronic  Corporation  Regulated  Power  Supply, 
Model  71,  was  used  as  a  probe  voltage  supply  when  a  non  time- 
varying  source  was  desired. 

A  Donner  Variable  Base  Function  Generator,  Model  3750,  was 
also  tried  in  an  attempt  to  automatically  make  plots  of  the 
natural  logarithm  of  probe  current  vs.  applied  probe  potential. 
Preliminary  attempts  were  unsuccessful  because  of  difficulty  in 
programming  a  suitable  logarithmic  function.   This  technique  was 
temporarily  abandoned  because  of  the  press  of  time.   It  does 
appear  quite  feasible,  nevertheless,  and  pursuit  along  this  line 
might  prove  to  be  quite  rewarding. 


* 

F.  W.  Moseley  Co.,  manufactures  a  logarithmic  attachment  for 
the  X  -  Y  recorder  which  enables  semi- logarithmic  plots  to  be 
made  without  need  of  a  function  generator.   They  also  manufacture 
a  curve  follower  attachment  to  the  X  -  Y  recorder  which  allows 
functions  to  be  drawn  on  paper  with  magnetic  ink;  the  pick-up 
coil  will  then  feed  out  the  function  desired  as  the  independent 
voltage  is  driven  onto  the  horizontal  axis.   A  system  such  as  ours, 
utilizing  this  technique,  would  require  two  recorders  --  one  to 
record,  one  to  generate  functions  --  both  being  driven  in  time 
synchronism. 
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l+.   Procedure . 

Prior  to  making  a  data  run,  the  vacuum  system  was  evacuated 
to  as  low  a  level  as  possible  --  in  the  range  of  10   mm  of  mercury. 
Then  the  system  was  baked  out  with  the  oven  at  4-50°  C  for  a  period 
of  between  twenty  four  and  thirty  six  hours.   During  this  time  the 
tubing  and  liquid  air  traps  external  to  the  oven  was  wrapped  in 
heating  tape  and  heated  to  350°  C.   Following  bake  out,  the  cathode 

was  heated  to  operating  temperature  with  the  filament  and  the  anode 

* 

was  outgassed  with  the  induction  heater.    The  system  was  then 

8       -9 
allowed  to  pump  down  to  between  5  x  10"'  and  10   mm  of  mercury. 

This  normally  took  24-  hours  after  completion  of  bake  out.   The 

system  was  then  filled  with  approximately  2  mm  of  spectroscopic- 

ally  pure  argon.   This  filling  procedure  was  repeated  several 

times  in  order  to  flush  the  tube  with  the  inert  gas.   With  argon 

Q 

in  the  system,  the  tube  could  be  evacuated  down  to  5  x  10   mm 
of  mercury  in  about  one  and  one  half  hours.   The  final  pressure 
was  obtained  by  discharging  the  tube,  measuring  the  pressure  with 
the  Langmuir  viscosity  gage,  and  correlating  it  with  the  plate 
voltage;  then,  evacuating  the  tube  in  a  small  increment  of  pressure, 
while  the  tube  was  still  discharging,  the  process  of  measuring 
pressure  and  correlating  with  tube  voltage  was  repeated  until 
pressure  -  tube  voltage  calibration  had  been  achieved.   The  final 
pressure  was  then  readily  obtained  by  evacuating  until  the  proper 
plate  voltage  was  obtained. 

This  bake  out  and  purification  procedure  is  tedious,  but  necess- 
ary.  For  example,  failure  to  be  careful  about  extreme  cleanliness 
results  in  a  pulsating  discharge  with  a  frequency  of  one  to  25  cps 
and  an  extremely  unstable  tube  characteristic  (i.e.,  the  tube  "runs 
away"  at  very  low  discharge  currents) . 
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An  arbitrary  interelectrode  separation  distance  of  about  M-2 
nun  was  selected  so  that  probe  measurements  could  be  made  with  the 
probe  at  increments  of  two  mm  along  the  axial  length  of  the  tube 
between  electrodes.   Thus,  about  twenty  different  positions  between 
anode  and  cathode  were  studied  at  approximately  equal  intervals. 

Electrode  positions  were  calculated  and  set  at  the  catheto- 
meter  as  follows:   The  cathetometer  was  set  at  its  reference  value 
for  the  cathode.   The  cathode  was  then  moved  using  the  magnetic 
coil  over  the  glass  surrounding  the  soft  iron  magnet  within  the 
tube  until  the  cathode  was  lined  up  with  the  reference  mark  on  the 
cathetometer.   The  cathetometer  was  then  traversed  to  the  probe 
position  where  its  value  was  recorded  (since  the  probe  position 
remained  fixed,  this  position  was  used  to  estimate  the  statistical 
accuracy  of  the  cathetometer  measuring  system) .  The  cathetometer 
was  then  traversed  to  the  anode  position  and  the  anode  moved  to 
coincide  with  the  reference  mark  as  described  for  the  cathode. 
All  cathetometer  readings  were  made  going  in  the  same  direction 
to  minimize  backlash  errors. 

A  tube  volt-ampere  characteristic  trace  was  then  made  and 
recorded  on  the  plate  circuit  X  -  Y  recorder  to  determine  the 
operating  range  of  tube  currents.   When  the  gas  in  the  tube  was 
very  clean  the  tube  remained  stable  up  to  200  milliamperes  dis- 
charge current.   Over  this  range,  then,  a  group  of  representative 
discharge  currents  were  selected  to  study  which  gave  suitable 
values  of  probe  current  (i.e.,  probe  currents  which  did  not  pro- 
duce voltages  in  excess  of  100  volts  to  the  operational  amplifiers 
of  the  analog  computer")  .   The  discharge  currents  selected  were 
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25,  50,  75,  100  and  125  milliamperes .   All  of  these  currents  were 
attainable  using  the  voltage  control  on  the  Kepco  Power  Supply, 
the  plate  circuit  resistance  remaining  fixed. 

Proceeding  next  to  the  probe  circuit,  calibration  procedures 
were  employed  to  minimize  all  electronic  errors.   The  amplifiers 
of  the  analog  computer  were  balanced  to  minimize  drift.   The 
function  multiplier  was  balanced  for  drift  and  accurate  suppression 
of  cross  products.   A  potentiometer  on  the  problem  board  was  used 
to  calibrate  the  maximum  peak-to-peak  excursion  of  the  triangle 
probe  voltage  supply  to  an  absolute  magnitude  of  150  volts  as 
traced  on  the  horizontal  axis  of  the  X  -  Y  recorder.   This  voltage 
was  then  biased  with  an  initial  condition  voltage  source  (set  with 
a  Donner  Null  VTVM)  to  give  the  maximum  positive  and  negative  ex- 
cursions desired.   Those  used  were  +20  volts  to  -130  volts  with  re- 
spect to  the  anode . 

Calibration  of  the  vertical  axis  voltage  corresponding  to  the 
probe  current  was  accomplished  with  a  problem  board  potentiometer 
after  the  signal  had  been  amplified  in  the  two  nominal  gain-ten 
stages.   This  potentiometer  was  then  adjusted  to  give  the  same 
magnitude  at  penpoint  on  the  X  -  Y  recorder  as  the  ammeter  reading 
probe  current.   This  had  to  be  accomplished  dynamically,  but  could 
be  done  to  quite  good  precision  (as  well  as  one  could  read  the 
meter)  at  very  slow  sweep  speed.   This  was  commonly  done  at  the 
maximum  ammeter  and  pen  excursion  when  both  momentarily  came  to 
rest.   Inertial  effects  were  considered  to  be  negligible. 
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Once  these  calibration  procedures  had  been  made,  data  runs 
were  started.   For  each  position  of  the  probe  a  linear  probe 
characteristic  curve  was  made  for  each  of  the  five  tube  discharge 
currents  cited  previously.   Blowup  traces  of  these  curves  were 
also  made  of  the  ion  (negative)  portions  of  the  trace.   The  paper 
was  then  changed  in  the  X  -  Y  recorder  and  the  current  squared  vs 
applied  probe  potential  plots  were  made  for  each  of  the  discharge 
currents.   Different  vertical  scales  were  used  for  electron  and 
ion  regions  of  the  curves  corresponding  to  those  used  for  the 
linear  traces.   Simultaneous  with  these  measurements,  the  plate 
circuit  volt-ampere  characteristic  was  being  monitored  and  traced 
giving  a  record  of  its  changes  during  the  probe  characteristic  run. 

Prior  to  each  new  position  of  the  probe  the  calibration 
procedures  cited  above  were  employed,  except  that  analog  computer 
amplifiers  were  only  balanced  approximately  hourly. 

A  run  normally  took  twenty  to  twenty  five  minutes  to  make  -- 
about  five  minutes  to  trace  the  curves,  and  about  fifteen  to  twenty 
minutes  to  reposition  electrodes  and  calibrate  equipment.   Twenty 
such  positions  between  anode  and  cathode  were  studied  using  this 
procedure  during  a  twelve  hour  period  yielding  four  hundred 
separate  probe  characteristic  traces. 

Reproducibility  of  curves  is  quite  good  within  the  limit  of 
electronic  errors  which  will  be  commented  upon  later  in  this  report. 
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A  subsequent,  more  limited,  data  run  was  made  for  ten  equi- 
spaced  probe  positions  (electrode  separation  40  mm,  argon  pressure 
6.9  microns  of  mercury,  discharge  current  50  ma),  for  the  purpose 
of  obtaining  second  derivatives  of  probe  current  with  respect  to 
probe  voltage.   Two  simple  RC  differentiating  circuits  in  cascade 
were  used  to  do  this  on  the  probe  current  signal,  (see  Fig.  6)  . 

Thse  derivative  curves  were  used  then  to  obtain  electron 
energy  distribution  functions. 

On  this  data  run,  ion  oscillations  were  measured  by 
oscilloscopic  photographs  and  electronic  counter  measurements. 
Electron  oscillations  were  measured  with  the  cavity  wavemeter 
coupled  to  the  dc  micro- ammeter. 

During  this  run  great  care  was  taken  to  evaluate  the  effect 
of  probe  contamination  by  retracing  a  curve  after  various  time 
intervals  and  noting  trace  deviations  from  the  original;  then 
after  cleaning  the  probe  by  ohmic  heating,  the  curve  was  again 
retraced  to  note  the  tendency  to  return  toward  the  original 
trace.   The  results  of  these  tests  were  quite  conclusive  and 
are  stated  later  in  Section  5,  Results. 
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5.   Results. 

Plots  of  two  representative  probe  positions  from  the  first 
run  (consisting  of  twenty  probe  positions  and  five  discharge 
currents)  are  shown  in  Figs .  7  through  13  for  the  following 

functional  relationships: 

2 

i  vs  V  ,  i    vs  V  ,  In  i  vs  V  ,  1^  vs  V^,  where  the  sub- 
script p  denotes  the  probe  parameter,  subscript  +  denotes  the  ion 
region  of  the  probe  characteristic  curve,  subscript  -  denotes  the 
electron  region  of  the  probe  characteristic  curve,  and  subscript 
d  denotes  the  tube  discharge  (or  plate)  parameter. 

A  summary  of  three  important  plasma  parameters  for  all 
twenty  probe  positions  between  anode  and  cathode  are  shown  in 
plots  as  indicated  below: 

(a)  n_  (electron  density)  vs  X  (cathode  -   probe  distance) 

2 
obtained  by  using  the  slope  of  i_   vs  V   curves  (Langmuir  Theory) 

is  shown  in  Fig.  14- . 

(b)  V    (space  potential)  vs  X,  obtained  from  In  ip  vs  V 
plots  at  the  first  positive  voltage  deviation  from  linearity  is 
shown  in  Fig.  15 . 

(c)  n._  vs  J  (discharge  current  density)  for  ten  probe 
positions  spaced  approximately  M-  mm  apart  to  correspond  with  the 
positions  used  during  the  second  data  run  is  shown  in  Fig.  16. 
These  data  are  perhaps  capable  of  being  correlated  with  the  re- 
sults of  the  second  data  run,  which  was  made  for  only  one  discharge 
current  (50  ma) ,  despite  the  fact  that  the  two  runs  were  made  at 
slightly  different  pressures. 
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A  second  data  run  was  made  at  50  ma  discharge  current  for 
10  probe  positions  spaced  at  4  mm  intervals  with  an  anode  -  cathode 
separation  distance  of  40  mm.   The  gas  pressure  was  6.9  microns  of 

Hg.   Results  of  this  run  are  indicated  below: 

2      2 

(a)  Plots  of  ip  vs  Vp,    di  /dV   ,  d  ip/dV  are  shown  in 

Figs.  17  a,  b   through  18  a,  b  for  two  representative  probe  posi- 
tions . 

(b)  Electron  energy  distribution  functions  for  all  ten  probe 
positions  are  shown  in  Figs.  19  through  28.   These  functions  were 
obtained  from  the  second  derivative  of  probe  current  with  respect 
to  probe  voltage  according  to  the  relation 

N0(eV)  =   2    pmVp  \  %      d2ip 

*e_ '       "dV7- 
p 

Where  N0(eV)  is  the  number  of  electrons  per  electron  volt  of  energy 
per  cubic  centimeter,  A  is  the  area  of  the  probe,  m  is  the  mass  of 
the  electron,  e  is  the  electronic  charge.   For  the  geometry  of  the 
probe  used  in  this  work,  the  function  reduces  to 

"1  ~i  i- 

N0(eV)  =  1.722  x  10  '  V2   d  in        _i   -3 

P      V  eV   cm 


2 

dV^ 

P 

Two  things  are  to  be  noted  concerning  these  derivatives  and 
associated  energy  distribution  functions.   First,  the  differentiat- 
ing circuit  was  not  perfect  --  in  particular  it  appears  that  the 
capacitor  was  not  unloading  quite  fast  enough  to  give  true  values 
when  the  derivative  was  decreasing.   This  error  was  relatively 
small  and  was  arbitrarily  adjusted  to  make  it  fit  the  first 
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derivative.   Second,  the  X  -  Y  plotter  developed  some  mechanical 
troubles  during  this  run  and  gave  a  magnitude  for  the  second 
derivative  of  approximately  1/3  of  the  true  value.   Approximate 
calibration  of  the  derivatives  was  achieved  by  using  another 
machine  which  was  reading  true  and  also  numerically  differentiating 
(well  realizing  the  errors  inherent  in  this  process)  the  i  vs  V 

XT  rr 

curve  to  establish  the  scale  factor.   The  second  derivatives  and 
the  distribution  functions  shown  in  Figs.  19  through  28  are  shown 
in  their  original  scale  as  they  came  from  the  machine.   The 

correction  factor  is  applied  to  the  constant  term,  i.e.  KCorrection 

13  %   d2' 
x  1.7  22  x  10   V       p    _   This  correction  factor  may  still  be 


P 


dV2 


P 
in  error  by  perhaps  ten  percent;  however,  these  results  are  shown 

because  the  technique  appears  valid  and  the  information  the  curves 
yield  appears  to  be  a  valuable  diagnostic  tool.   The  actual  magni- 
tudes really  are  not  critical  for  diagnostic  purposes;  shapes  and 
changes  in  shape  in  going  from  one  probe  location  to  another  are 
important  because  these  factors  show  how  the  energy  of  the  particles 
is  being  transferred. 

(c)   A  plot  of  n_  vs  X  is  shown  in  Fig.  29.    These  values 
were  obtained  by  integrating  the  energy  distribution  functions,,   By 
this  method  the  density  could  be  divided  between  the  high  and  low 
energy  groups  and  changes  in  density  groupings  noted  as  the  probe 
position  was  changed.   The  integration  was  carried  out  numerically 
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using  the  trapezoidal  rule.    Results  were  checked  with  a  plan i 
meter  and  agree  within  three  per  cent.   Fig  i.   M7  shows  a  typical 
distribution  function,  and  the  technique  employed  to  integrate  it 
and  obtain  the  various  densities. 

(d)   A  plot  of  n_  vs  X  is  shown  in  Fig.  30.   These  data  were 

obtained  by  the  conventional  Langmuir  method  of  using  the  slope  of 

2 
i   vs  V  plots.   Comparison  of  Fig.  29  and  Fig.  30  shows  that  both 

plots  have  the  same  general  shapes,  and  essentially  the  same  peaks 
and  gradients  are  obtained  by  both  methods.   The  absolute  magni- 
tudes vary  by  as  much  as  30%.   This  may  be  attributable  to  the 
scale  factor  assumed  for  the  second  derivative  method;  but,  possibly 
part  of  the  difference  (perhaps  ten  per  cent)  may  exist   in  comput- 
ing the  slopes  by  the  Langmuir  method  since  correction  terms  were 
not  used,  and  only  a  relatively  few  points  were  used  to  determine 
these  slopes.   Since  the  main  purpose  in  obtaining  these  values 
was  to  check  the  approximate  magnitude  of  those  values  obtained  from 
the  distribution  functions,  great  significance  is  not  attached  to 
these  small  differences.   Even  though  absolute  magnitudes  of  the 
densities  computed  from  the  distribution  functions  are  not  certain3 
the  magnitudes  relative  to  each  other  are  thought  to  be  quite  good 
and  this  fact  is  used  later  to  correlate  density  with  the  predicted 
Langmuir  frequency  and  to  make  a  statistical  estimate  of  the 

The  trapezoidal  rule  gives  an  approximate  value  to  the  area  under 
an  irregular  curve,  i.e.,  Area  2=    ^(y^  +  y2  +  „ „ .  + yn     +  yn/2)  . 
In  our  integration  we  used  a  ^  of  two  electron  volts  giving  a  fair 
approximation  to  the  integral.   Using  smaller  increments  gives  an 
increasingly  better  result.   Simpson's  rule  is  more  accurate.   We  used 
the  trapezoidal  rule  because  it  is  suitable  for  making  a  running  tabu- 
lar integral.  Simpson's  rule  is  not  conducive  to  this  treatment. 
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standard  deviation  for  the  values  of  density  obtained. 

(e)  A  plot  of  relative  electron  density  vs  X  is  shown  in 

Fig.  31.   These  data  were  obtained  from  In  i  vs  V  curves  by 
&  P      P         * 

using  the  current  at  space  potential  and  the  electron  temperature, 

(f)  A  plot  of  relative  electron  density  of  high  and  inter- 
mediate energy  beam  electrons  is  shown  in  Fig.  32.   These  data 
were  obtained  by  Waymouth's  technique  of  analysis  of  i  vs  V 
probe  characteristic  curves  to  obtain  the  energy  of  the  group 
Medicus'  technique,  based  upon  a  diffusion  law  which  is  dependent 
upon  pressure,  was  used  for  relative  magnitudes.   Our  pressure 
range  was  significantly  different  (Lower)  from  that  of  Medicus* 
development;  hence,  the  results  are  good  only  to  an  order  of 
magnitude . 

(g)  A  plot  of  eV  (average  electron  energy)  vs  X  is  shown  in 
Fig.  33.   These  data  were  obtained  from  the  average  energy  of  the 
two  electron  groups.   The  average  energy  of  each  group  being 
obtained  from  the  relation 

eV? 

(  v3/2    ah 

eV  =    J       _   d(eV) 

eV,  c 

L  dV 

P 


f2    V       d2- 


P 


d(eV) 


These  values  are  shown  on  the  distribution  functions  in  Figs.  19 
through  28. 
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(h)  A  plot  of  T  (electron  temperature)  vs  X  is  shown 
in  Fig.  3S.  These  data  were  obtained  from  the  slope  of  the 
linear  portions  of  the  In  i  vs  V  curves  after  the  current 
attributable  to  the  high  energy  electrons  had  been  subtracted  out. 

(i)   A  plot  of  Te  vs  X  for  both  high  and  low  energy  electron 

groups  is  shown  in  Fig.  36.   These  data  were  obtained  from  the 

average  energy  of  the  two  electron  groups  by  using  the  relation 

eV  =  3/2  kT   .   This  is  strictly  true  only  if  the  distribution  is 

Maxwellian.   The  distributions  obtained  were  close  to  this  for  the 

low  energy  groups  and  the  relation  was  used  for  the  high  energy 

groups  as  well  for  want  of  a  better  relation  to  use.   This  plot  is 

essentially  the  same  as  Fig.  33,  since  both  were  obtained  from 

the  same  data.   The  ordinates  differ  only  by  a  constant. 

(j)   A  plot  of  Vcri  vs  X  is  shown  in  Fig.  37.   These  data 
t>p 

were  obtained  by  taking  the  maximum  of  the  first  derivative  for 
positive  values  of  V  as  being  space  potential.   This  corresponds 
to  the  inflection  point  of  the  i  vs  V  curve,  or  the  first  devia- 
tion from  linearity  on  the  In  i_  vs  V  plot. 

P     P 

(k)   A  plot  of  V   vs  X  is  shown  in  Fig.  38.   These  data  were 
t>p 

obtained  by  taking  the  first  deviation  from  linearity  on  the  In  i„ 
vs  Vp  plot  as  being  space  potential.   Agreement  between  this  figure 
and  Fig.  37  is  quite  good.   More  faith  is  placed  in  the  values  ob- 
tained from  the  first  derivative  because  these  values  are  easier  to 
read  precisely  than  the  deviation  from  linearity  on  the  logarithmic 
plot. 

(1)   A  plot  of  E  (electric  field  strength)  vs  X  is  shown  in 
Fig.  39.   These  data  were  obtained  by  taking  the  changes  in  space 

potential  as  a  function  of  distance  (X)  from  Fig.  37. 
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Plasma  oscillations  in  two  frequency  ranges  werp  found. 
Kilomegacycle  signals  detected  are  believed  to  be  electron 
oscillations  and  kilocycle  signals  are  believed  to  be  ion 
oscillations . 

High  frequency  oscillations  of  47  2  Mcs  and  520  Mcs  were 
found  at  positions  20  mm  and  24-  mm,  respectively,  from  the  cathode 
for  the  6.9  micron,  50  ma  discharge.   Computing  the  density  by 
of  Langmuir's  theory  relating  the  density  to  the  slope  of  iz      vs 
V  curves,  we  obtain  the  plot  of  total  density  shown  in  Fig.  30a, 
below.   The  densities  corresponding  to  the  observed  frequency, 
calculated  by  using  Langmuir's  relation  f  =  8980  n2,  are.  shown 
as  dashed  lines  on  this  plot. 
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Fig.  30a.   Electron  density  as  a  function  of  cathode  -  probe 

distance.  Data  obtained  from  slopes  of  ±^   vs  Vn  plots 
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Neither  of  these  densities,  corresponding  to  the  observed 
frequencies,  correspond  with  measured  points  on  the  total  density 
plot,  but  instead,  represent  plasma  oscillation  frequencies  lower 
than  would  be  obtained  from  total  density  values. 

However,  if  we  make  use  of  the  density  plot  obtained  from  the 
energy  distribution  functions  and  also  use  Waymouth's  technique  to 
obtain  the  distribution  of  high  energy  beam  electrons,  we  obtain 
the  plot  of  densities  shown  in  Fig.  29a,  below.   By  using  these 
techniques  we  are  able  to  break  our  electron  concentrations  down 
into  groups  which  are:   high  energy  beam  electrons  (obtained  by 
Waymouth's  technique)  (see  p.  23,  Section  2,  Previous  Work), 
intermediate  energy  beam  electrons  (obtained  from  energy  distri- 
bution functions)  and  random  plasma  electrons  (obtained  from 
energy  distribution  functions) . 
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Fig.  29  a.  •   Total  electron  density  (from  distribution  functions) 
O  Plasma  Density  (from  distribution  functions) „ 
X  Intermediate  energy  beam  electrons  (from 
distribution  functions) 
—  High  energy  beam  electrons  (Waymouth's  technique) . 
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If  we  apply  the  relationship   cited  by  Kharchenko,  el   i 
(31)  (described  on  page  26  ,  Section  2,  Previous  Work) . 

CD         fo  =    fp 


1   - 


2  v  i 

s 


2 

v 


where  S  is  the  thermal  plasma  velocity  and  v   is  the  beam 

o 

electron  velocity  and  put  this  in  terms  of  energy,  we  have 
(2) 


'b 
where  eV  is  the  average  thermal  plasma  energy  and  eVV  is  the 

IT 

average  energy  of  the  beam  electrons. 

We  know  f  ,  eV  ,  and  eVV  experimentally.   Let  us  now  calculate 
the  equivalent  density  corresponding  to  f   and  see  where  this  fits 
on  the  density  plot,  Fig.  29a. 

Consider  first  f  =  520  Mcs. 


(3)      eV 


P 


eVb 


=  5.5/22  =0.25 


intermediate 


((+)   and,    f   =   (1  -  0.25)  2   520  Mcs 

IT 


f   =   451  Mcs 
P 


Realizing  this  was  developed  for  vQ  >^^>  S 
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The  equivalent  density  is 

/  6    \    2  q  3 

(5)  ,  =451x10  =      2.51  x  10     particles/cm 

Similarly,    for   fQ      =   472   Mcs, 

(6)  eV 

_  =      5/50      =      0.1 

bhigh 

h 

(7)  and,    f        =      (1   -    0.1)       472  Mcs 

f  =      447   Mcs. 

P 

The   equivalent  density   is 

(8)  n2_  =  2W?   x  1Q9  particies/cm3o 

Referring  to  Fig.  29a,  densities  n^  and  np   (and  their 
corresponding  frequencies)  are  indicated  by  dashed  lines.   It  is 
apparent  that  these  densities  correspond  with  the  curve  repre- 
senting the  plasma  density,  not  total  density . 

In  obtaining  the  data  for  Fig.  29a3we  were  beset  by  a 
dilemma.   We  have  previously  mentioned  that  the  densities  ob- 
tained by  the  second  derivative  method  were  in  error  by  some 
constant  factor  caused  by  improper  machine  calibration.   However 3 
we  are  confident  that  the  relative  values  of  these  densities  are 
rather  good  with  respect  to  each  other.   Now  the  absolute  values 
of  densities  obtained  by  the  Langmuir  method  are  probably  accurate 
to  about  10%.   In  order  to  use  our  density  data  obtained  by  the 
second  derivative  method  to  best  advantage  we  have  determined  a 
scale  factor  to  be  applied  to  these  densities  (obtained  by  the 
second  derivative  method)  by  taking  ratios  of  each  corresponding 
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ordinate  of  Fig.  29  (second  derivative  data)  and  Fig.  30 
(Langmuir  data) .   We  have  used  the  average  of  these  ratios 
to  correct  both  curves  to  the  same  absolute  magnitude. 

Then,  by  statistically  analyzing  the  ten  density  ratios  so 
obtained,  we  are  able  to  compute  a  number  which  gives  a  certain 
measure  of  confidence  in  the  values  of  density  plotted  in  Figs. 
29a  and  30a.   We  call  this  number  s  to  denote  a  standard  devia- 
tion  ratio,  and  plot  is  as  0""  =   +  s  n_  £   particles/  cm  on  Figs. 
29a  and  30a,  for  the  total  density  and  plasma  density  curves, 
s  was  calculated  to  bt  +  0.17.   n_  ^    is  used  to  denote  any  ordinate 
to  the  density  plots  and  a~  the  standard  deviation  in  particles/ 
cm  . 

Now  it  is  quite  interesting  to  note  that  neither  density 
corresponds  to  the  Langmuir   frequency  within  the  confidence 
limits  when  we  attempt  to  infer  it  from  the  total  density  plot. 
However,  when  we  explain  our  observed  frequencies  by  correlating 
them  with  the  plasma  density  curve,  we  note  that  these  densities, 
n^   and  np  ,  fall  within  one  standard  deviation  of  this  curve,  if 
we  consider  that  part  of  the  plasma  density  curve  at  the  point  of 
disintegration  of  the  beam  density  curve. 

The  conclusions  to  be  drawn  from  this  analysis  are  brought 
out  in  the  next  section  of  this  report. 

It  is  worthwhile  to  note  that  only  the  472  Mcs  oscillation  had 
a  large  amplitude  (2  ma  of  detector  current) .   The  520  Mcs  oscilla- 
tion was  much  smaller  in  amplitude  (0.01  to  0.02  milli-micro  amp 
of  detector  current)  .   We  associate  the  strong  M-7  2  Mcs  oscillation 
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with  the  intermediate  energy  beam  electrons  (eVu  =  22  electron  voll^O 
and  the  520  Mcs  oscillation  with  the  high  energy  beam  electrons 
(eVb   =  50  electron  volts) .   Note  also  that  the  density  of  this 
group  of  high  energy  beam  electrons  is  about  one  order  of  magnitude 
smaller  than  either  the  plasma  electrons  or  the  intermediate  energy 
beam  electrons . 

Other  high  intensity  electron  oscillations  were  detected  at 
higher  currents  much  more  readily  than  at  50  ma  discharge  current. 
Unfortunately,  this  was  discovered  so  late  as  to  preclude  taking 
correlating  probe  data.   The  results  of  these  measurements,  for 
one  probe  position,  are  shown  in  tabular  form  below. 


Discharge   Current 
(ma) 

Observed 
Frequency 
(Mcs) 

Approximate 
Computed 
Frequency 
(Mcs) 

Inter 
Sign; 
(Det€ 
Curre 

3    x 

isity  of 
il 

;ctor 

?nt   -    amp) 

100 

703 

710 

lO"11 

105 

704.5 

720 

1.2 

_Q 

x  10   ' 

110 

708.5 

730 

4  x 

io-7 

115 

71. 

737 

4    X 

nf7 

120 

713.3 

746 

4    X 

IO"7 

125 

717.5 

752 

6   x 

io"7 

For  this  probe  position  no  oscillation  signals  were  detected 
below  100  ma  discharge  current.   This  suggests  that  a  threshold 
effect  for  interaction  may  be  present  as  a  function  of  current 
density. 

* 

By  interpolating  in  Fig.  14  using  a  pressure  correction  factor, 
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Fig.  M-0  is  a  plot  of  current  density  (Jd)  vs  observed  fn 

2 
quency  squared  (f  )  showing  the  linear  relation  of  these  values. 

The  other  type  of  oscillation  found  was  a  low  frequency  sig- 
nal of  about  100  Kcs ,   The  ion  frequency  predicted  by  calculation 
from  the  density  of  the  intermediate  energy  beam  electrons  is 
135  Kcs.   Hence,  these  oscillations  are  believed  to  be  ionic 
in  nature.   They  existed  throughout  the  length  of  the  tube  (for 
the  50  ma  discharge  at  6.9  microns  pressure),  maintaining 
essentially  the  same  frequency;  however,  the  intensity  of  the 
signal  did  vary  with  position  in  the  tube.   The  results  of  these 
intensity  variations  are  shown  in  Fig.  3M-,  where  I  (intensity  in 
volts)  is  plotted  vs  X.   It  is  also  significant  to  note  that  the 
oscilloscopic  traces  of  these  oscillations  show  a  more  fixed 
adherence  to  a  steady  fundamental  frequency  and  waveform  on  the 
cathode  side  of  the  discharge  where  the  high  and  intermediate 
energy  electron  beams,  coming  from  the  cathode,  were  in  evidence. 
On  the  anode  side  of  the  discharge  (in  the  region  after  which 
scattering  has  occurred)  these  oscillations  show  a  wider  spectrum 
of  frequencies  (varying  between  90  and  110  Kcs)  and  the  intensities 
of  signals  and  waveshapes  are  qui Le  variable  when  sampled  in  time. 
Fig.  M-l  is  an  oscillographic  photograph  of  an  oscillation  pattern 
on  the  cathode  side  of  the  meniscus  region.   Fig.  M-2  is  an  oscillo- 
graphic photograph  of  an  oscillation  pattern  obtained  on  the 
anode  side  of  the  ,ianiscus  region. 
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Visual  and  photographic  observation  of  the  discharge  \ 
made  in  effort  to  try  to  discover  a  meniscus,  a  scattering  region 
and/or  a  diverging  beam  of  electrons.   No  meniscus  was  observed 
at  low  currents;  but,  visual  observation  seemed  to  indicate  a 
region  of  maximum  disturbance  in  the  beam,  and  at  about  the  same 
location,  a  tendency  for  the  beam  to  diverge  from  the  shape  dis- 
played in  emerging  from  the  cathode.   The  region  in  which  this 
phenomenon  appeared  is  indicated  on  Figs.  29  through  39  as  the 
cross-hatched  region.   Figs.  43  through  45  are  photographs  and 
sketches  of  the  discharge  which  show  the  beam,  its  tendency  to 
diverge,  and  apparent  meniscus  formation.   Note  that  meniscus 
formation  appears  to  grow  with  increasing  current. 

In  attempting  to  correlate  frequency  measurements  made  at 
higher  currents  during  the  second  data  run  (50  ma  discharge  at 
6.9  microns  of  pressure)  with  the  current-density  data  of  the 
first  run  (five  different  currents  at  5  9  microns  of  pressure), 
an  interesting  fact  (overlooked  as  being  trivial  when  it  occurred) 
was  noted.   Referring  to  Fig.  14-  for  a  discharge  current  of  125  ma, 
at  a  position  ten  mm  from  the  cathode,  a  peak  density  occurs.   The 
minimum  of  this  gradient  occurs  at  12  mm  from  the  cathode.   Referring 
to  the  original  probe  characteristic  curve  corresponding  to  this 
point  (Fig.  7,  125  ma  discharge  current),  it  was  noted  that  consid- 
erable irregularity  of  the  ion  portion  of  this  curve  was  observed 
in  a  probe  voltage  range  of  -40  to  -100  volts  with  respect  to  the 
anode.   The  phenomenon  manifested  itself  in  a  step  increase  in 
probe  current  of  about  0.004  ma  which  corresponded  to  a  step 
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increase  of  about  0  .  <4-  volts.   This  distortion  in  probe  chara<  i  er- 
istic lasted  out  to  about  -  100  volts  at  which  point  i1  t<  nded  to 
return  to  its  normally  expected  shape,  but  with  some  noticeable 
change  in  slope.   At  the  time  of  occurrence,  the  peculiarity  was 
noticed  and  the  curve  retraced  several  times,  each  time  giving  a 
similar  result,  although  the  exact  voltage  at  which  the  step  in- 
crease in  current  occurred  did  not  coincide  exactly,  it  held  to 
within  a  20  volt  region.   Sweep  speed  of  probe  voltage  was  in- 
creased with  no  change  in  effect  --  in  fact  the  speed  of  sweep 
was  increased  to  about  50  cps  (vs  the  normal  0.05  cps)  and  the 
trace  photographed  on  the  oscilloscope.   The  same  sort  of 
phenomenon  showed  thereon  except  that  this  trace  seemed  to  show 
oscillating  signals  riding  on  the  probe  characteristic  from  about 
-20  volts  out  to  the  maximum  excursion  of  -130  volts.   Further 
examination  of  the  electron  portion  of  the  probe  characteristic 
curve  disclosed  that  this  was  the  first  probe  characteristic  curve 
(for  this  current) ,  in  moving  the  probe  from  cathode  to  anode,  in 
which  a  section  having  a  constant  slope  in  the  positive  current 
region  (denoting  a  region  separating  two  energy  groups,  ref. 

discussion  of  Waymouth's  work,  pg.  23   ,  Section  2,  Previous 

* 
Work)  did  not  exist. 

A  similar  occurrence  took  place  for  the  same  discharge  current 

at  10  mm  from  the  cathode.   This  disturbance  to  the  probe  was  less 

pronounced,  amounting  to  a  change  in  slope  on  the  linear  i   vs  V 

curve  and  a  distinct  deviation  from  linearity  on  the  ±        vs  VL  curve 

P+     P 

See  Fig.  M-6  which  illustrates  this.  cf.  Fig.  8,  in  which  only 
one  electron  group  is  present. 
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Still  considering  this  125  ma  discharge,  it  was  noted  at 
the  time  of  taking  the  data  that  the  position  10  mn>      the 
cathode  was  just  barely  on  the  cathode  side  of  the  diverging 
beam  region  (observed  visually)  of  the  discharge.   The  electron 
portion  of  this  probe  characteristic  curve  showed  a  constant 
slope  region  denoting  two  distinct  groups  of  electrons  still 
present  at  this  position  in  the  tube.   The  position  12  mm  from 
the  cathode  was  inside,  i.e.,  on  the  anode  side,  of  the  diverg- 
ing beam.   The  meniscus  is  believed  to  have  been  present  about 
10  mm  from  the  cathode. 

Similar  occurrences  were  evident  on  the  probe  character- 
istic curves  at  14  mm  from  the  cathode  for  the  100  ma  discharge  and 
at  22  ram  (see  Fig. 8)  from  h^  cathode  for  the  75  ma  and  125  ma,  dis- 
charges at  probe  voltages  of  about  -100  volts  with  respect  to 
the  anode .   These  points  corresponded  with  a  peak  density 
gradient  in  each  case  (see  Fig.  14)  . 

In  light  of  what  was  discovered  on  the  second  data  run 
(50  ma  discharge  at  6.9  microns  pressure)  it  now  appears  clear 
what  happened  in  these  cases  during  the  first  data  run.   These 
disturbances  to  the  probe  characteristic  are  associated  with  the 
same  region  in  which  (a)  electron  oscillations  were  detected., 
(see  Fig.  29) .   (b)   ion  oscillation  intensity  reached  a  maximum 
and  ion  oscillation  waveshape  became  most  distorted,  (see  Fig. 
(c)  beam  electron  concentration  largely  disappeared  and  merged  into 
one  random  distribution  (see  Figs.  23  and  2M-)  ,  and  (d)  either  a 
maximum  or  a  minimum  density  gradient  was  present  (see  Fig.  14)  . 
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The  changes  in  slope  as  far  negative  as  -100 
believed  to  be  caused  by  electrons  acquiring  enough 
the  interacting  region  in  the  vicinity  of  the  meniscus  (wh< 
oscillations  are  believed  to  originate)  so  thdt  they  can  ov<  rcoiro 
this  amount  of  retarding  field  and  become  collected  by  the  probi 
There  are,  it  appears,  distributions  of  electrons  in  these  inter- 
acting regions  having  maximum  energies  as  high  as  100+  electron 
volts.   The  step  increase  in  current  shown  on  rig,,         .ieved 
to  be  a  monoenergetic  group  of  electrons  having  an  energy  of 
about  50  electron  volts  with  a  variability  of  about  +  10  electron 
volts  when  sampled  in  time. 

A  distinction  can  be  made  with  respect  to  the  effect  on  the 
probe  characteristic  as  it  is  affected  by  the  density  gradient 
A  large  disturbance,  i.e.,  a  step  increase  in  probe  current,  and 
a  small  disturbance,  i.e.,  a  change  in  slope^was  noted  where  the 
position  was  associated  with  a  minimum  of  the  density  gradient 
(corresponding  with  Emeleus'  findings  that  oscillations  reach  a 
maximum  intensity  at  the  bottom  of  a  density  gradient)  .      Only  a 
small  disturbance,  i.e.,  a  change  in  slope  of  the  probe  character- 
istic curve,  was  noted  where  the  position  corresponded  to  a  raa> 
mum  of  a  density  gradient  (in  agreement  with  Emeleus'  findings 
that  oscillations  are  initiated  at  or  near  peaks  in  the  density 
gradient) . 
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The  effect  of  probe  contamination  was  studied  to  det- 
whether  cathode  sputtering  was  causing  the  probe  to  coli«  d  Li 
pure  material  and  thus  change  its  collecting  propertii  ■-„   A 
curve  was  traced  with  a  freshly  cleaned  probe  (by  ohm;-  hi  i  I  ng 
to  a  cherry-red  condition)  then  allowed  to  sit  immediately  in 
front  of  the  cathode,  perhaps  2  mm  away,  (the  worse  possible 
position) ,  for  a  period  of  15  minutes  to  one  ha  I  houi  md  th< 
the  same  characteristic  curve  was  retraced,   A  slight  change  in 
the  characteristic  was  noticeable;  however,  when  the  interval  of 
time  was  reduced  to  the  same  time  as  was  used  for  making  a  set  of 
data  curves,  the  effect  of  probe  contamination  on  the  character- 
istic curves  was  entirely  negligible.   This  is  probably  attribut- 
able to  the  very  low  evaporation  rate  of  material  from  the  cath( 

2 
given  as  less  than  0.8  micrograms  of  barium  per  cm  per  hour  by 

the  manufacturer. 

Miscellaneous  tests  were  conducted  to  determine  the  effects 
of  electronic  errors  associated  with  the  measurement  of  probe 
characteristic  curves.   Probably  the  best  estimate  of  the  rel ; 
bility  of  the  probe  characteristic  curves  was  obtained  by  allow- 
ing the  curve  to  be  retraced  15  to  20  times  to  note  the  drift 
tendency  due  to  the  amplifiers  and  associated  electronic  circuitry. 
Results  of  these  measurements  showed  that  the  maximum  deviation  to 
be  expected  was  of  the  order  of  +  0.002  ma  of  current  for  electron 
currents  of  the  order  of  o '1  ma.   Curves  were  also  re-run  after  an 
elapsed  time  of  days  and  no  significant  change  in  characteristic 
was  noticed  except  when  the  gas  was  obviously  dirty  because  of 
long  periods  of  operation  (up  to  several  weeks) . 
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6 .   Conclusions . 

The  results  of  our  second  data  run  (50  ma  discharge  at  6 
microns  pressure  in  argon)  as  shown  in  Fig,  29   tend  to  substanti- 
ate that  plasma  electron  oscillations  originate  at  the  peak  of  a 
density  gradient  and  grow  in  amplitude  down  a  density  gradient  as 
originally  proposed  by  Allis  (30)  and  observed  experimentally  by 
Emeleus  and  co-workers  (33) . 

Electron  oscillations  were  detected  only  at  positions  in  the 
tube  immediately  after  divergence  of  the  electron  beam  had  occurred. 
These  oscillations  could  be  associated  with  a  peak  density  gradient 9 
and  the  observed  frequency  of  these  electron  oscillations  bore  a  re- 
lationship to  the  predicted  Langmuir  frequency  when  correlated  with 
the  plasma  density  (vs   the  usual  total  density) .   The  reason  for 
this  will  become  clear  later  in  discussing  a  possible  mechanism 
for  oscillation  excitation. 

Ion  oscillations  were  detected  throughout  the  length  of  the 
tube  with  a  fundamental  frequency  of  about  100  Kcs.   Fig.  34-  shows 
that  the  intensity  of  these  oscillations  fell  to  a  minimum  just 
prior  to  the  point  where  the  total  density  gradient  reached  a 
maximum  (see  Figs   29  and  34)  and  then  the  oscillation  intensity 
rose  rapidly  to  a  maximum  where  the  total  density  gradient  was  a 
maximum.   This  region  is  presumed  to  be  that  where  particle  inter- 
action is  greatest  and  where  energy  exchange  is  taking  place  and 
thus  exciting  electron  oscillations .   This  region  is  outlined  by 
cross-hatching  on  Figs.  29  rhrough  39.   The  fact  that  the  photo- 
graphic traces  of  the  ion  oscillations  across  this  region;,  in 
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going  from  cathode  to  anode,  show  considerable  regularity  on  th< 
cathode  side  (see  Fig.  41)  and  more  intensity  and  more  irregularity 
in  wave  shape  and  frequency  on  the  anode  side  (si    Lg.  L(-2) 
indicates  that  in  this  region,  strong  microfields  are  present 
which  are  disturbing  the  ion  oscillations.   The  change  in  intensity 
of  the  ion  oscillations  indicates  that  ions  are  taking  part  in  the 
energy  exchange  mechanism  by  losing  energy  as  they  move  out  of  the 
turbulent  region  toward  the  cathode . 

Disturbances  to  the  probe  characteristic  curves  noted  on  the 
first  data  run  being  associated  with  maximum  or  minimum  density 
gradients  points  strongly  to  the  conclusion  that  strong  fields 
and  oscillations  were  present,  and  considerable  energy  exchange 
was  taking  place,  at  these  positions  in  the  tube,   Disturbances 
out  to  -100  volts  suggests  that  some  electrons  were  acquiring 
better  than  100  electron  volts  of  energy   These  fast  electrons 
must  be  generated  in  the  interacting  region  since  the  maximum 
energy  that  could  be  gained  from  the  discharge  tube  would  be 
of  the  order  of  the  tube  voltage  drop,  about  20  electron  volt- 

Maximum  disturbances  were  noted  at  positions  corresponding 
to  a  minimum  of  the  density  gradient  in  agreement  with  Alii 
contention  that  oscillations  grow  in  intensity  down  a  decreasing 
density  gradient.   Smaller  disturbances  to  the  probe  character- 
istic curve  were  noted  at  positions  corresponding  to  a  maximum 
of  the  density  gradient  in  accord  with  Allis'  idea  that  plasma 
oscillations  originate  at,  or  near,  a  peak  density  gradient.   At 
higher  discharge  currents,  multiple  density  gradients,  (i.e.,  two 
peaks)  were  present  (see  Fig.  1M-) .   Probe  disturbances  were 
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associated  with  each  of  these.   This  fact  sugges'   that  a  more 

complex  interaction  is  taking  place  with 

of  two  such  interacting  regions,  one  following  ti     u?r. 

Examination  of  the  distribution  f unc I      shown  in  Figs.  19 
through  28  show  that  two  groups  of  electrons  are  presenl  on  the 
cathode  side  of  the  interacting  (meniscus)  region .   On  the  am 
side  of  the  interacting  region,  only  one  group  of  electrons  is 
present.   Further,  by  employing  Waymouth's  analysis  of  the  probe 
characteristic  curve,  a  third  higher  energy  group  of  electrons  is 
also  found  to  be  present,  but  these  electrons  are  either  much  less 
numerous  or  of  too  high  an  energy  to  show  on  the  energy  distribu- 
tion function.   We  have  earlier  denoted  these  higher  energy  groups 
as  an  intermediate  energy  group  of  beam  electrons  and  a  high  energy 
group  of  beam  electrons,  respectively. 

These  beam  electrons  emerge  from  the  region  in  front  of  the 
cathode  and  interact  with  the  plasma  electrons  at  a  region  in  the 
tube  where  visual  divergence  of  the  beam  and  the  aforementioned 
disturbances  in  the  probe  characteristic  curves  occur. 

The  evidence  of  both  data  runs  has  suggested  to  us  that  the 
mechanism  exciting  the  oscillations  is  associated  with  these  beam 
electrons,  because  all  of  these  associated  phenomena  are  indicative 
of  a  region  of  instability  in  which  local  fluctuating  microfields 
of  high  strength  would  be  likely  to  exist. 

In  the  light  of  this  evidence,  the  authors  believe  that  the 
beam  electrons,  upon  reaching  the  plasma,  are  strongly  interacting 
with  the  plasma  (when  beam  current  densities  are  adequate) 9  pro- 
ducing local  instabilities  which  excite  plasma  oscillations. 
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Ion  oscillations   may  well  play  a  role   in  the   excitation  process, 
or  may  only  be    a   by-product   of   the   electron  beam  excitat  i 
mechanism. 

The  theory  of  such  a  mechanism  has  been  developed  by  Akhiezer 
and  Fainberg  and  is  discussed  by  Kharchenko,  Fainberg,  Nikolayev, 
et  al.,  (31).   The  important  points  in  their  discussion  are  sum- 
marized in  Section  2,  Previous  Work,  page  26. 

As  previously  mentioned,  there  are  two  possible  types  of 
instabilities,  convective  and  alsolute„   In  our  discharge,  con- 
vective  instability  is  believed  to  have  occurred  since  oscillations 
were  detected  only  in  a  very  small  region  associated  with  the 
initiation  of  the  instability,  and  none  were  detected  on  the  anode 
side  of  the  point  of  initial  detection  where  one  would  expect  th< 
oscillation  intensity  to  increase  as  the  oscillation  moved  down  a 
density  gradient.   Possibly  the  spread  in  velocities  as  the  beam 
electrons  become  randomized  in  their  initial  interaction  with  the 
plasma  tends  to  remove  the  instability  as  suggested  by  Kharchenko, 
et  al.,  (31).   This  may  be  the  reason  that  oscillations  are  not 
detected  other  than  in  the  region  of  initial  scattering. 

Perhaps  as  a  bit  of  speculation,  we  might  mention  that  the 
possibility  of  modified  "runaway"  electrons  (31,  1,  36)  being 
generated  and  initiating  the  instability  should  not  be  entirely 
ignored.   Although  "runaway"  electrons  have  not  been  observed 
(to  our  knowledge)  in  discharge  tubes  such  as  ours,  conditions 
suitable  for  their  production  appear  to  be  met  in  our  tube  and 
some  of  the  phenomena  observed  might  well  be  attributed  to  their 
presence.   For  example,  the  sudden  termination  of  the  instability 
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and  damping  out  of  oscillations  are  what  one  would  expect  if  I 
excitation  was  produced  as  a  result  of  "runaway"  electrons. 

Or  possibly,  electrons  traversing  the  intei  icting  '  gion  r 
be  given  a  pulse  of  energy  by  the  oscillations         ge  as 
"runaway"  electrons.   These  "runaway:-"  would  not  interact  with 
other  particles  because  their  collision  pro! 

less  at  their  greater  energy,  thus,  they  would  completely  emerge 
from  the  plasma  and  finally  give  up  some  of  the '   exi   -  energy 
by  collision  with  walls  of  the  tube. 

One  other  point  is  worthy  of  discussion  at  this  time.   On  the 
second  data  run,  oscillations  were  detected  at  higher  dischargf 
currents  having  signals  of  considerably  greater-  amplitude  th 
those  observed  at  50  ma  discharge  current.   All  efforts  at  detect 
oscillations  at  discharge  currents  lower  than.  50  ma  were  unsuccess- 
ful.  Examination  of  the  density  plot  from  the  first  data  r 
(Fig.  14)  shows  that  no  sharp  gradients  exist  for  the  25  and  50 
ma  discharges.   Sharp  gradients  do  exist  for  75  ma  and  above.   In 
the  light  of  this  observation  and  the  experience  obtained  on  the 
second  data  run,  it  is  believed  that  current  density  (in  order  to 
produce  a  sufficiently  strong  interacting  beam) ,  and  pressure  are 
important  in  establishing  the  instability  and  hence  oscillatic 

In  particular,  had  we  operated  at  slightly  higher  pressures 
and/or  slightly  higher  currents  ,  it  is  believed  thai   ;<  illations 
would  have  been  found  more  readily  and  that  a  closer  approach  to 
an  ideal  plasma  (i.e.,  having  a  Maxwellian  distribution  of  energie 
would  have  been  obtained  since  a  smaller  mean  I    path  would  have 
produced  more  randomization  of  electron  energit- 
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The  possibility  also  exists  that  we  were  operating  in  a 
current  density  -  pressure  region  which  was  just  adequate  to  allow 
a  beam  to  form  and  move  some  distance  before  interacting  with 
plasma.   Operation  at  higher  current  densities  has  shown  u 
the  meniscus  region  moves  toward  the  cathode  as  current  is 
Further  experimentation  should  reveal  the  optimum  current  density 
pressure  relationship  to  use  in  attempting  to  isolate  the        i< 
mechanism. 

The  results  of  this  investigation  are  based  upon  only  limited 
data,  some  of  which  is  questionable  in  magnitude 

per  cent  in  error  in  ordinates  of  the  energy  distribution  functions) 
because  of  the  equipment  troubles  at  the  end  of  our  work.   •   • 
all  data  were  qualitatively  good,  at  least.  Much  more  data  nee 
to  be  collected  before  definite,  positive  conclusions  may  be  drawn. 
Nevertheless,  the  correlation  between  all  of  the  evidence  presented 
is  quite  convincing  and  seems  to  substantiate  the  tentative  con- 
clusions which  have  been  drawn. 


The  area  of  the  cathode  used  in  this  work  is  ff    cm  ;  hf 
dividing  any  of  our  discharge  currents  by  "IT  will  give  the  < 
density . 
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7 .   Recommendations . 

In  keeping  with  our  findings  and  observations,  the  following 
recommendations  for  further  work  are  ma'' 

(a)  Continue  probe  studies  and  oscillation  detection  at 
higher  pressures  and/or  higher  currents,  perhaps 
varying  pressure  as  a  parameter  to  determine  changes 
in  the  form  of  distribution  functions. 

(b)  Establish  the  pressure-current  density  range  in  which 
oscillations  are  detected  most  readily,  then  attempt 
to  find  a  simple  oscillation-density  gradient  rela- 
tionship which  will  present  the  least  complicated 
structure  initially. 

(c)  Once  suitable  oscillations  are  found,  make  probe 
studies  and  frequency  measurements  of  both  ion  and 
electron  oscillations  to  determine  where  in  the  tube 
they  arise  and  attempt  to  correlate  these  measure- 
ments with  visual  observation,  energy  distribution 
functions,  and  perhaps  also  some  spectroscopic 
observations . 

(d)  Study  the  region  in  which  oscillations  appear  to 
arise  in  great  detail;  i.e.,  make  probe  and  oscilla- 
tion measurements  around  this  region  utilizing  very 
small  increments  of  movement  of  the  probe,  say  0.5 
mm  or  less,  in  an  effort  to  isolate  the  mechanism 
of  excitation. 
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Figure    1 
General    Experimental    Arrangement 
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Figure  2 
Characteristic  Amplifying  and  Recordii 
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Figure  3 
Power  Supplies  and  Probe  Circuitry 
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Figure  Ua 
Discharge  Tube 
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DISCHARGE    TUBE 


Figure  kb 
Discharge  Tube 
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Figure  29.   Electron  density  as  a  function  of  Cathode  ■=  Probe  position, 

Obtained  from  energy  distribution  functions, 
•  --  Total  density 
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Figure  30.   Total  electron  density  as  a.function  of  Cathode 
position.   Obtained  from  the  slope  of  i    vs.   V   curves. 
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Figure  31.   Relative  total  electron  density  as  a  function  of  Cathode  - 
Probe  position.   Obtained  from  In  i   vs.   V  curves  near  space  potential, 
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Figure  32.   Relative  electron  density  of  beam  electrons  as  a  function 
of  Cathode  -  Probe  distance. 

0  --   Intermediate  energy  group 

X. --  High  energy  group 
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Figure  33.   Average  energy  of  electrons  in  electron  volts  as  a  function 
of  Cathode  -  Probe  distance.   Obtained  from  energy  distril  tion  functions 
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Figure  34.   Ion  oscillation  intensity  as  a  function  of  Cathode  •=  Probe 
distance. 
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Figure  35.   Electron  temperature  in  degrees  Kelvin  as  a  function  of 

Cathode  -  Probe  distance.   Obtained  from  In  i   vs.   V   curves. 
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Figure  36.   Average  electron  temperature  in  degrees  Kelvin  as  a  function 
of  Cathode  -  Probe  distance.   Obtained  from  energy  distribution  functions 
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Figure  37.   Space  potential  in  volts  as  a  function  of, Cathode  - 

distance.   Obtained  from  maximum  of  first  derivative  of  i   vs. 
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Figure  39.   Electric  field  strength  as  a  function  of  Cathode 
distance.   Obtained  from  data  of  Figure  37. 
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Figure  40.   Current  density  as  a  function  of  the  square  of  the  observed 
frequency  of  electron  oscillation  for  the  probe  located  immediately 
after  the  region  of  suspected  maximum  interaction. 
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Fig  41a.   Several  oscilloscope  traces 
of  ion  oscillations. 


Fig  41b.   Oscilloscope  trace  of  ion 
oscillation  when  sampled  in  time. 


Fig  41c.   Oscilloscope  trace  of  ion 
oscillation  when  sampled  in  time. 


Fig  41.   Oscilloscope  traces  of  ion  oscillations  on  the 
cathode  side  of  the  meniscus  region.   Discharge  current 
50  ma.   Pressure  6.9  microns  of  mercury.   Cathode-anode 
snacing  40  ran.   Probe  position  12  mm  from  cathode  - 
oscilloscope  sweep  speed  5y«sec/cm.   Amplitude  1  volt/cm. 
■ote  smaller  amplitude,  more  fixed  waveform  and  more  steady 
frequency  compared  with  Fig  42. 
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Fig  42a.   Several  oscilloscope  traces 
of  ion  oscillations. 


Fig  42b.   Oscilloscope  trace  when 
sampled  in  time. 


Fig  42c.   Oscilloscope  trace  when 
sampled  in  time. 


Fig  42.   Oscilloscope  traces  of  ion  oscillations 
On  the  anode  side  of  meniscus  region.   Cathode-anode 
spacing  40  mm.   Probe  position  20  mm  from  cathode. 
Discharge  current  50  ma.  Pressure  6.9  microns  of 
mercury.   Oscilloscope  sweep  speed  5/^sec/cm. 
Amplitude  1  volt/cm.   Note  the  greater  amplitude, 
distortion  of  wave  form  and  spreading  of  fundamental 
frequency  compared  with  Fig  41. 
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Figure  hj    a.   Photograph  of  50  ma  discharge 
in  argon  at  6.9  microns  of  mercury. 
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Figure  k3    b.   Sketch  of  50  ma  discharge  in 
argon  at  6.9  microns  of  mercury. 
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Figure  kU-   a.  Photograph  of  120  ma  discharge 
in  argon  at  6.9  microns  of  mercury 
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Figure  kk   b.   Sketch  of  120  ma  discharge  in 
argon  at  6.9  microns  of  mercury. 
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Figure  k5    a.   Photograph  of  greater  than  150  ma 
discharge  in  argon  at  6.9  microns  of  mercury. 


Anode 


Cathode 


Visible  meniscus 


X    (mm) 

Figure    1+5   b.       Sketch   of  greater   than   150   ma 
discharge    in   argon   at    6.9   microns    of   mercury. 


Figure   k5 
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Fig   h7&.    Energy   Distribution  Function.  Argon 
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It  is  convenient  to  express  this  computation  in  tabular  form 
in  which  a  running  integral  may  be  set  up  and  the  integral 
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V  N0(eV)  d(eV)  computed  at  the  same  time  for  the  purpose 
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of  evaluating  the  average  energy,  eV,  of  a  group  of  electrons, 


Figure  U7.   Example  of  electron  density  computation  and  average 
energy  computation  from  energy  distribution  function. 
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These  computations  are  shown  in  a  tabulation  on  the 
next  page  in  Figure  47b. 


Figure  k7 .   Example  of  electron  density  computation  and  average 
energy  computation  from  energy  distribution  function. 
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